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Abstract  of  Dissertation  Presented  to  the  Graduate  School  of 
the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

PREDICTIVE  MODELING  OF  HIGH-CURRENT  OUTPUT  RESISTANCE  AND 
THERMAL  EFFECTS  IN  BIPOLAR  JUNCTION  TRANSISTORS 

BY 

SANG-GUG  LEE 
August  1992 

Chairman:  Dr.  R.  M.  Fox 

Major  Department:  Electrical  Engineering 

This  dissertation  concerns  predictive  modeling  of  high- 
current  output  resistance  and  thermal  effects  in  bipolar 
junction  transistor  (BJT) . Modern  BJTs  tend  to  operate  with 
saturated  carrier  velocity  in  the  collector  space-charge 
region  (SCR) . A physical  analysis  of  the  effects  of  velocity 
saturation  on  the  output  resistance  of  BJTs  is  presented, 
especially  for  high  current  levels.  Physical  analyses  show 
that  when  the  collector  current  density  approaches  a critical 
value,  the  Early  voltage  can  increase  significantly  due  to 
carrier-velocity  saturation  as  the  resulting  base  push-out 
partially  offsets  base-width  modulation.  This  effect  is  also 
demonstrated  in  simulations  with  PISCES  and  MMSPICE.  The 
predicted  peak  in  Early  voltage  cannot  be  observed  in  low- 
frequency  measurements  due  to  the  thermal  effects.  The 
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effects  of  self-heating  on  BJT  behavior  are  demonstrated 

S 

through  measurement  and  simulation.  Most  affected  are  the 
small-signal  parameters  Y22  and  Y12  • A frequency-domain 
solution  to  the  heat-flow  equation  is  presented  which  applies 
to  any  rectangular  emitter  geometry.  The  model,  although 
simple  enough  for  CAD,  predicts  thermal  spreading  impedance 
with  good  accuracy  for  a wide  range  of  frequencies.  A 
practical  method  is  presented  for  extracting  thermal 
spreading  resistance  of  BJTs.  The  method  uses  the  output 
resistance  or  reverse  transconductance  as  the  temperature- 
sensitive  parameter.  Measurement  can  be  made  in  the  time  or 
frequency  domain.  General  insights  into  local-heating 
effects  on  the  large-  and  small-signal  behavior  of  integrated 
circuits  are  provided,  and  the  circuit  conditions  under  which 
self-heating  causes  significant  errors  in  simulation  are 
identified.  A thermally  improved  BJT  model  is  used  as  a 
verification  tool  in  studying  the  implication  of  thermal 
effects  on  integrated  circuits. 
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CHAPTER  1 
INTRODUCTION 


Analog  circuits  such  as  operational  amplifiers  typically 
employ  active  loads  to  achieve  high  gains.  The  voltage  gain 
and  frequency  response  of  such  circuits  thus  depends  on  the 
transistors'  small-signal  output  resistance.  In  most 
contemporary  bipolar  junction  transistors,  high  electric 
fields  occur  in  the  collector  junction  space-charge  region 
due  to  the  high  epi-layer  doping  or  high  reverse  bias 
voltage.  These  high  fields  tend  to  saturate  the  free-electron 
drift  velocity,  and  thereby  affect  output  resistance. 

In  Chapter  2 of  this  dissertation,  a study  is  done  of 
the  effects  of  carrier-velocity  saturation,  especially  for 
high  current  levels,  on  the  output  resistance  characteristics 
of  BJTs . Physical  and  numerical  analyses  show  that  when  the 
collector  current  density  approaches  a critical  value,  the 
Early  voltage  can  increase  significantly  due  to  carrier- 
velocity  saturation  as  the  resulting  base  push-out  partially 
offsets  base-width  modulation.  This  effect  is  demonstrated  in 
simulations  with  PISCES  and  MMSPICE.  Measurements  suggest 
that  thermal  feedback  tends  to  mask  this  effect  at  low 
frequencies.  The  discrepancies  between  analysis  and  low- 
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frequency  measurements  led  us  to  study  the  thermal  effects  in 
bipolar  junction  transistors. 

Existing  electronic  circuit  simulation  programs,  both 
empirical  (Gummel-Poon,  SPICE  2 and  3)  and  physics-based 
(PISCES,  MEDUSA,  BIPOLE,  MMSPICE),  generally  model  the 
integrated  circuit  as  if  all  the  devices  within  the  circuit 
are  at  a constant  temperature.  In  fact,  the  power  dissipated 
by  the  elements  of  the  integrated  circuit  can  raise  the 
instantaneous  temperature  of  the  element  itself  and  other 
elements  above  ambient  temperature.  Increases  in  device 
temperature  can  lead  to  significant  increases  in  terminal 
currents  and/or  decreases  in  base-emitter  voltage  depending 
on  the  bias  conditions.  The  thermally  induced  variations  can 
also  lead  to  significant  changes  in  y -parameters , especially 
Y12  and  Y22- 

There  has  been  a considerable  amount  of  work  done  on  the 
modeling  of  thermal  feedback  in  BJTs  [Sp58,  St59,  Mu64,  Jo70, 
Me77,  Da88,  De89] , yet  these  effects  are  neglected  in  most 
existing  device  and  circuit  simulators.  One  probable  reason 
for  this  is  the  empirical  approach  which  is  popular  in  BJT 
modeling.  For  example,  the  Early  voltage,  commonly  used  to 
model  output  resistance,  actually  drops  at  high  currents,  yet 
a constant  value  is  still  commonly  used.  A second  reason  is 
that  thermal  effects  are  considered  to  be  a problem  unique  to 
power  transistors. 
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The  temperature  rise  in  a device  is  proportional  to  the 
thermal  resistance,  defined  as  the  total  temperature  rise 
caused  by  a unit  power  step.  A packaged  transistor  can  be 
characterized  by  two  thermal  impedance  components,  viz.,  the 
thermal  spreading  impedance  controlling  the  fast  temperature 
variation  and  the  package  thermal  impedance  controlling  the 
slow  process.  The  package  thermal  impedance  is  dependent  on 
the  thermal  properties  of  the  packaging  material  and  cause 
coupling  of  all  the  devices  on  the  chip.  The  thermal 
spreading  impedance  is  concerned  with  the  geometry  of  the 
transistor.  In  the  past,  the  contribution  of  the  thermal 
spreading  resistance  to  the  total  thermal  impedance  has  been 
negligible.  'However,  circuit  and  device  integration  levels 

V 

continue  to  increase  with  time.  The  resulting  smaller  device 
size  leads  to  higher  thermal  spreading  impedance  so  that  it 

tends  to  dominate  the  total  thermal  impedance.  At  the  same 

' s 

time,  the  demand  for  high  performance  has  led  to  higher  power 
densities  in  these  devices.  The  result  of  these  trends  is 
towards  larger  temperature  rises  in  device.  Therefore,  in 
modern  high-performance  circuits,  the  device  thermal 

resistances  and  power  densities  are  sufficient  to  result  in 

v 

substantial  temperature  rises  within  a circuit.  Also,  new 
isolation  technologies  (e.g.,  trench  isolation,  dielectric 
isolation,  and  SOI)  tend  to  make  thermal  effects  more 
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significant  due  to  the  poor  thermal  conductivity  of  silicon 
dioxide . 

; A well-known  attribute  of  the  flow  of  heat  from  a small 
source  into  large  media  is  its  effective  confinement  to  a 
region  less  than  a few  times  the  lateral  dimension  of  the 
source  [Ho58] . For  transistors  dissipating  large  power, 
heating  within  the  elements  of  the  integrated  circuit  can 
cause  considerable  chip  temperature  gradients.  In  this  case, 
the  power  dissipation  not  only  influences  the  performance  of 
the  device  itself  but  leads  to  thermal  coupling  between 
devices.  This  type  of  thermal  coupling  is  well  known  and  in 
most  cases  they  can  be  avoided  by  using  symmetrical  layout, 
placing  in  close  proximity,  or  arranging  critical  pairs  of 
transistors  as  cross-connected  quads.  For  modest  powers,  even 
adjacent  devices  are  thermally  decoupled  from  each  other 
[Fo91a]  so  that  localized  heating  of  the  BJT  is  intrinsic  to 
the  device  itself  and  fundamental.  This  dissertation  focuses 
on  this  intrinsic  self-heating  effect  which  cannot  be  avoided 
through  changing  the  layout  of  the  IC. 

Joy  and  Schlig  [Jo70]  produced  an  excellent  time-domain 
model  for  the  thermal  response.  In  their  time-domain 
derivation,  Joy  and  Schlig  [Jo70]  integrated  a point-source 
heat-flow  equation  over  the  volume  of  the  collector  space- 
charge  region  (SCR)  and  its  image.  The  result  is  a closed- 
form  expression  for  the  impulse  response  at  any  point  in  the 
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semiconductor.  The  time  integral  of  this  expression  gives  the 

response  to  a unit  step  in  the  collector  power.  The  integral 

is  complicated  and  does  not  have  a closed-form  solution.  The 
step  response  at  t = «*>  gives  the  thermal  resistance  Rth- 

Their  model  is  too  CPU-intensive  for  use  in  circuit 
simulation.  Furthermore,  thermal  feedback  tends  to  cause  more 
significant  errors  in  the  small-signal  behavior  of  a circuit 
than  in  the  transient  response.  A simple,  yet  scalable  and 
predictive  frequency-domain  model  of  thermal  spreading 
impedance  was  needed  for  small-signal  thermal  modeling. 

Recently  several  workers  [Vo89,  Ye91]  have  demonstrated 
ways  to  model  the  thermal  feedback  in  circuit  simulation 
assuming  thermal  model  parameters  are  available.  While  the 
thermal  spreading  resistance  is  the  dominant  part  of  the 
total  thermal  impedance  in  scaled-down  devices, 
characterizing  these  devices  is  very  difficult.  Based  on  Joy 
and  Schlig's  model,  thermal  time  constants  of  the  scaled 
devices  are  on  the  order  of  nanoseconds.  The  standard 
electrical  method  for  thermal  impedance  measurement  is  the 
“emitter-only  switching"  method,  and  it  is  adapted  for 
characterizing  power  transistors  [Oe76] . In  this  technique,  a 
large  heating  current  is  first  applied  to  establish  thermal 
equilibrium.  This  current  is  then  turned  off  and  the 
temperature  is  inferred  from  the  resulting  transient  response 
of  Vbe»  measured  at  a small  fixed  emitter  current.  The  initial 
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response  is  controlled  by  both  electrical  and  thermal 
spreading-impedance  effects.  This  makes  it  difficult  to 
determine  the  original  temperature,  and  prevents  accurate 
extraction  of  the  thermal  spreading  resistance  Rth-  Therefore, 
there  was  a need  for  more  reliable  characterization  method 
for  thermal  parameters . 

It  was  shown  [De89]  that  neglecting  intrinsic  local 
heating  effects  caused  errors  up  to  4%  in  the  predicted  delay 
of  a simple  BJT  logic  gate.  Local  heating  can  severely 
degrade  the  performance  of  certain  types  of  analog  integrated 
circuits.  The  types  of  circuits  that  are  subject  to  thermal 
effects  are  those  circuits  which  require  high  precision  such 
as  D/A  converters,  instrumentation  amplifiers,  analog 
multipliers,  logarithmic  amplifiers,  current  mirrors,  and 
precision  voltage  references,  and  those  types  of  circuits 
which  experience  large  amounts  of  power  dissipation  such  as 
voltage  regulators,  audio  and  servo  amplifiers.  Local  heating 

V— 

can  also  cause  errors  in  op-amp  practice  as  it  can  lower  the 
low-frequency  open-loop  gain,  cause  distortion,  change  the 
phase  response  as  well  as  the  offset  voltage  and  common-mode- 
rejection  of  differential  amplifiers.  Therefore,  it  is 
important  for  designers  to  have  good  insight  into  the  circuit 
implications  of  local-heating  effects,  as  even  well-executed 
designs  can  fall  prey  to  these  effects  [To90] . Knowledge  of 
the  implications  of  thermal  effects  in  integrated  circuits 
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not  only  provides  the  design  principles  to  reduce  thermal 
effects  in  IC  design  but  can  be  used  to  exploit  them  for 
better  IC  design. 

This  dissertation,  then,  is  concerned  with  predictive 
modeling  of  thermal  effects  in  bipolar  junction  transistors 
and  the  circuit  implications  of  thermal  effects.  This  work 
will  provide  a good  insight  on  the  high  current  output- 
resistance  characteristics  and  a general  understanding  and 
prediction  of  thermal  effects  in  BJT  circuits.  The  main 
contributions  made  in  this  work  are  as  follows: 

(1)  predictive  modeling  of  the  effects  of  carrier-velocity 
saturation,  especially  for  high  current  levels,  on  the 
output-resistance  characteristics  of  BJTs; 

(2)  a simple  and  yet  predictive  frequency-domain  thermal 
impedance  model; 

(3)  measurement  techniques  for  characterizing  thermal 
impedance ; 

(4)  general  insights  into  local-heating  effects  on  the  DC 
and  small-signal  performance  of  integrated  circuits  and 
identif cation  of  the  circuit  conditions  under  which 
self -heating  causes  significant  errors  in  simulation. 
Chapter  2 presents  a physical  analysis  of  the  effects  of 

velocity  saturation  on  the  output  resistance  of  BJTs. 
Analytical  expressions  for  normalized  BJT  output  resistance 
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are  derived  for  an  idealized  structure.  MMSPICE  and  PISCES-II 
simulations  are  also  provided. 

In  Chapter  3,  the  effects  of  self -heating  on  BJTs  are 
demonstrated  through  measurement  and  simulation.  This  chapter 
will  show  that  local-heating  causes  much  larger  errors  in  BJT 
small-signal  parameters  than  large-signal  parameters.  A 
thermal  impedance  model  is  derived  that  predicts  thermal 
corrections  to  BJT  small-signal  parameters  based  on  a 
solution  of  the  heat-flow  equation  in  the  frequency  domain. 
The  model  is  applicable  in  principle  to  any  emitter  geometry. 
Polynomial  approximations  provide  simplified  expressions  for 
rectangular  emitters.  These  expressions  are  simple  enough  for 
compact  modeling  in  SPICE. 

In  Chapter  4,  practical  methods  are  presented  for 
extracting  the  thermal  spreading  resistance  of  BJTs 
separately  from  package  thermal  impedance.  The  methods  use 
the  output  resistance  or  the  reverse  transconductance  as  the 
temperature-sensitive  electrical  parameter.  Measurement 
techniques  in  both  time  and  frequency  domain  are  introduced. 

The  thermal  impedance  model  described  in  Chapter  3 has 

N 

been  implemented  in  SPICE  DC  and  AC  analysis  by  D.  Zweidinger 
[Zw92] . In  Chapter  5,  the  improved  SPICE  BJT  model  [Fo91c]  is 

l 

introduced  which  includes  a brief  description  of  how  the 
thermal  models  are  incorporated  in  the  DC  and  AC  routines, 
and  discusses  the  three  levels  of  thermal  modeling  available 
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to  the  user.  In  addition,  a series  of  simulations  are 
compared  to  measurements  to  demonstrate  the  improved  accuracy 
afforded  by  this  thermally  enhanced  BJT  model.  DC  thermal 
effects  are  considered  in  two  aspects,  magnification  of  the 
increments  in  collector  current  due  to  the  mismatch  in 
collector-base  voltage  and  collector-current  deviation  from 
near-exact  exponential  dependence  on  base-emitter  voltage.  On 
the  small-signal  side,  the  relations  between  the  thermal 
effects  on  the  input  and  output  resistances  and  the  series 
emitter  (Re) , base  (Rb) , and  collector  (Rc)  resistances  are 
discussed.  The  implication  of  discrepancies  in  the  input  and 
output  resistance  to  other  small-signal  aspects  are  also 
discussed.  Simulation  results  of  a few  circuit 
configurations,  based  on  modified  and  ordinary  SPICE,  are 
provided  for  verification  and  demonstration  purposes.  Design 
strategies  to  avoid  significant  thermal  problems  are  also 
discussed  in  this  chapter. 

x. 

In  Chapter  6,  the  main  accomplishments  of  the 
dissertation  are  summarized,  and  suggestions  of  areas  for 


future  work  are  made . 


CHAPTER  2 

THE  EFFECTS  OF  CARRIER -VELOCITY  SATURATION  ON  HIGH-CURRENT 

BJT  OUTPUT  RESISTANCE 

2.1.  Introduction 

Analog  circuits  such  as  operational  amplifiers  typically 
employ  active  loads  to  achieve  high  gains . The  voltage  gain 
and  frequency  response  of  such  circuits  thus  depend  on  the 
transistors'  small-signal  output  resistance.  In  most 
contemporary  bipolar  junction  transistors,  high  electric 
fields  occur  in  the  collector  junction  space-charge  region 
due  to  the  high  epi-layer  doping  or  high  reverse  bias 
voltage.  These  high  fields  tend  to  saturate  the  free-electron 
drift  velocity,  and  thereby  affect  output  resistance. 

This  chapter  presents  a study  of  the  effects  of  carrier- 
velocity  saturation,  especially  for  high  current  levels,  on 
the  output  resistance  characteristics  of  BJTs.  Analytical 
expressions  for  normalized  BJT  output  resistance  are  derived 
for  an  idealized  structure  as  an  aid  to  qualitative 
understanding.  MMSPICE  [Je89b]  and  PISCES-II  [Pi85] 
simulations  are  also  provided. 
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2.2.  Analysis 

In  a recent  letter,  Roulston  [Ro90]  presented  an 
improved  model  for  the  Early  voltage,  Va,  in  NPN  transistors 

with  very  narrow  base-widths  (less  than  0.1|im).  Roulston 
concluded  that,  due  to  carrier-velocity  saturation,  a 
substantial  improvement  in  Early  voltage  is  expected  in 
narrow-base  transistors  compared  to  estimates  based  only  on 
base  doping  level  and  width.  Liou's  [Li90]  correction  to  this 
work  led  to  the  following: 

Va  = Vao  [ 1 + Dn  F^/Vg  Wbo  1 (2.1) 

where  Vao  = Wbo/ (dWBo/dVcE)  is  the  conventional  expression  for 
Early  voltage  (neglecting  velocity  saturation) , Dn  is  the 
electron  diffusion  constant,  Fc  is  an  empirical  factor  to 
account  for  built-in  fields  and  bandgap  narrowing  in  the  base 
region,  vs  is  the  electron  saturation  velocity,  and  Wbo  is  the 
quasi-neutral  base  width.  For  a silicon  bipolar  transistor 
with  a base  width  of  0.04  |im,  Dn  = 17  cm2/V-s,  and  Fc  = 1, 
this  model  predicts  Va  = 1.4  Vao  • 

With  carrier- veloci ty  saturation,  the  carrier 
concentration  in  the  space-charge  region  (SCR)  is 
proportional  to  the  collector  current,  and  the  positions  of 
the  SCR  boundaries  become  functions  of  the  carrier 
concentration.  Roulston' s model  accounts  only  for  a reduction 
in  the  concentration  gradient  due  to  the  nonzero  carrier 
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concentration  at  the  SCR  boundary.  From  (2.1),  this  leads  to 
significant  increase  in  Va  only  for  very  narrow  base  widths. 
As  is  well  known,  at  the  level  of  collector  current  where  the 
SCR  boundaries  become  strong  functions  of  Jc,  base-width 
modulation  due  to  increments  in  Vcb  tends  to  be  offset  by  base 
push-out  (Kirk  effect)  [Ki62]  as  the  resulting  increase  in  nc 
(and  Jc)  redistributes  the  electric  fields  in  the  SCR.  This 
offsetting  mechanism,  which  was  neglected  in  [Ro90] , can  lead 
to  a much  more  significant  increase  in  small-signal  output 
resistance  at  high  currents,  even  for  ordinary  base-widths. 

Following  is  a derivation  of  this  increase  in  output 
resistance  as  represented  by  the  Early  voltage.  The  Early 
voltage  [Ge78]  can  be  considered  as  the  output  resistance 
normalized  by  the  collector  current:  Va  = rQIc  _ vce  where  rQ 

is  the  output  resistance  and  Ic  is  the  collector  current.  To 
simplify  the  mathematics,  base  and  buried-layer  dopings  are 
assumed  constant,  and  the  buried-layer  doping  is  assumed  much 
greater  than  the  base  doping,  which  is  much  greater  than  the 
epi-doping.  The  base  terminal  is  assumed  to  be  driven  by  a 
constant  voltage  source  (Fig.  2.1)  and  the  emitter  junction 
SCR  width  is  assumed  to  be  independent  of  the  collector 
current.  The  base-collector  bias  voltage  and  the  resulting 
electric  field  are  assumed  to  be  high  enough  to  saturate  the 
electron  velocity  even  near  the  SCR  edge.  By  this  assumption 
we  are  excluding  the  possibility  of  ohmic  quasi-saturation 
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Fig.  2.1.  One-dimensional  NPN  transistor  structure. 
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even  though  in  practice  the  transistor  might  enter  that 
operating  region  at  low  Vcb  • Ohmic  quasi-saturation  is 
excluded  from  this  analysis  since  it  can  be  treated  without 
considering  velocity  saturation. 

Three  operating  conditions  are  considered:  (Case  1) 
forward-active  with  part  of  the  epitaxial  layer  quasi- 
neutral, (Case  2)  forward-active  with  the  epi-layer  entirely 
space-charged,  and  (Case  3)  non-ohmic  quasi-saturation 
[Je89a] . 


2.2.1  Case  1 

Figure  2.2a  shows  the  electron  concentration  and 
electric  field  distribution  for  the  case  of  forward-active 
operation  with  part  of  the  epi-layer  quasi-neutral  (QN) . In 
the  SCR,  the  current  flows  predominantly  by  drift.  Assuming 
velocity  saturation,  the  collector  current  density  Jc  can  be 

expressed  as 

Jc  = q ncvs  (2.2) 

where  nc  is  the  electron  concentration  in  the  space-charge 
region  required  to  support  the  current.  Equation  (2.2)  and 
the  equality  of  charges  on  both  sides  of  the  collector 
junction  lead  to  the  relation 


xB  (qNA  + Jc/vs)  = xc  (qNEPi  - Jc/vs) 


(2.3) 
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Fig.  2.2.  Typical  electron  concentration  and  electric  field 

distribution  for  (a)  forward-active  operation  with 
part  of  the  epitaxial  layer  quasi-neutral,  (b) 
forward  active  operation  with  the  epi-layer 
entirely  space-charged,  Jc  < Jepi.  (c)  forward 
active  operation  with  the  epi-layer  entirely  space- 
charged,  Jc  > Jepi.  (d)  non-ohmic  quasi-saturation 
operation. 
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Fig.  2.2  continued 
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where  xb  is  the  penetration  of  the  SCR  into  the  base,  xc  is 
the  penetration  of  the  SCR  into  the  epi-layer,  Nepi  is  the 
epi-layer  doping  concentration,  and  Na  is  the  base  doping 
concentration.  Poisson's  equation  can  be  solved  in  the 
collector  junction  SCR  to  yield  the  following  expression: 


Vcb  + <|>c 


qNA  + 


+ 


|qNEPi  - 


(2.4) 


where  Vcb  is  the  collector  junction  reverse  bias,  <j>c  is  the 
collector  junction  built-in  potential,  and  es  is  the 

semiconductor  permittivity.  In  (2.4)  the  potential  drop 
across  the  quasi-neutral  portion  of  the  epi-layer  is 
neglected  for  simplicity.  The  effect  of  this  approximation  on 
output  resistance  will  be  discussed  later.  From  (2.3)  and 
(2.4),  xb  can  be  written  in  terms  of  Vcb  and  Jc  as: 


xB 


/\l  2£sNepi  ( Vcb  + foe ) / 1 - Jc/JEpi  _ /~1~-  Jc/ Jepi 

V qNA  (Nepi  + Na)  v 1 + J q/ Ja  V 1 + Jc/ Ja 


where 


xB0 


V2  £gNEpi  (Vcb  + ^c) 

qNA  (NEpi  +Na) 


(2.6) 


Jepi  - qNgpiVg,  and  Ja  = QNavs.  xbo  is  the  SCR  penetration  into 
the  base  which  would  result  from  neglecting  velocity 
saturation  (i.e.,  using  the  depletion  approximation). 

The  quasi-neutral  base-width  (Wb)  with  carrier-velocity 

saturation  can  be  expressed  by 
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WB  = Wjc  - Xb 


(2.7) 


while  the  base  width  neglecting  velocity  saturation  (Wbo)  is 
expressed  by 

Wbo  = Wjc  - xbo  • (2.8) 

Figure  2.3  shows  a plot  of  xb/xbo  as  a function  of  collector 
current  density  Jq.  As  can  be  seen  from  Fig.  2.3,  as  Jc 
approaches  Jepi,  xb  decreases  rapidly.  Therefore,  dWB/dVcE 
differs  significantly  from  the  value  dWBo/dVcE  found  in 
Roulston  [Eqn.  2.1]. 

Jc  can  be  expressed  as 

Jc  = KqDn  [n (0 ) -nc]/WB  (2.9) 

where  K is  an  empirical  factor  that  allows  for  drift  current 

and  n(0)  is  the  electron  concentration  at  the  edge  of 
emitter-base  SCR  in  the  base  region  (Fig.  2.2a).  Using  (2.2), 
(2.9)  and  the  definition  of  the  Early  voltage,  Va  can  be 

expressed  as 

VA  = - [WB  + KDn/vs]  (dVCE/dWB)  . (2.10) 

Therefore,  from  (2.5),  (2.7),  and  (2.10),  we  obtain 


0.5  xB 


KDn  + 
vsWB 


Esvs  (Vcb  + ^c)Jc 
xbWb(Jc  + Ja)2  .. 


(2.11) 
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Fig.  2.3.  Plot  of  the  ratio  xb/xbo  as  a function  of  collector 

current  density  Jc  for  Na  = lxlO17  cm-3  and  Nepi  = 
lxlO15  cm-3. 
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Similarly,  an  expression  for  the  Early  voltage  neglecting 
velocity  saturation  can  be  found  from  (2.8),  (2.9),  and 

(2.10),  assuming  nc=0,  as 

— , _ WBo  ( Vcb  + <|>c ) 

V AO  — 

0 . 5 Xbo  . ( 2 . 12 ) 

Equation  (2.11)  is  similar  in  form  to  (2.1)  except  for 
the  addition  of  a third  term  in  the  brackets.  The  last  two 
terms  have  only  a small  effect  on  Va  except  for  very  narrow 
base-widths.  For  Jc  > 0.1  JEpi,  however,  (2.11)  predicts  more 
significant  increases  in  Va  due  to  the  ( 1-Jq/ Jepi ) 1 1 2 
dependence  of  xb  in  the  factor  Wb  (Vqb+^c)  / (0 . 5 xb)  in  (2.11). 

Equation  (2.11)  is  valid  only  up  to  the  value  of  Jc  for 
which  the  base-collector  SCR  extends  all  the  way  to  the 
buried  layer  for  a given  value  of  Vcb#  i.e.,  when  xc  = Wepi. 
The  upper  limit  of  validity  for  (2.11)  can  be  found  from 
(2.3)  and  (2.5)  with  xc  replaced  by  Wepi: 


dc  Ixc=Wepi 


J EPI 

+ Na  I Xbq  \ 2 
NEPi  \WEpx/ 


1 - 

XC0  \ 2 

WepiI  - 

(2.13) 


Xco 

where 


/ 2esNA  (VCb  + <1>c) 
V qNEPI  (NEpx  + NA) 


(2.14) 


is  the  value  of  xc  obtained  from  the  depletion  approximation. 
For  typical  BJTs,  (Na/Nepi)  (xbo/Wepi)2  « 1 • Thus,  when  Wepi  >> 
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xco ' (2.11)  is  valid  for  values  of  Jc  very  close  to  Jepi/  so 
that  the  effect  of  carrier-velocity  saturation  on  Va  can  be 
very  large. 


2.2.2.  Case  2 

As  Jc  increases  above  JclXc  = wEPI  the  collector-base  SCR 

begins  to  uncover  charge  in  the  buried  layer.  The  resulting 
electron  concentration  and  electric  field  distribution  are 
shown  in  Fig.  2.2b.  As  the  current  increases  above  Jepi,  the 
electron  concentration  in  the  epi  rises  above  Nepi,  leading  to 
field  and  electron  concentration  distributions  as  shown  in 
Fig.  2.2c.  Using  charge  neutrality  and  solving  Poisson's 
equation  over  the  collector  junction  SCR,  xb  for  Case  2 

(Figs.  2.2b  and  2.2c)  can  be  written  as 


From  (2.7),  (2.10),  and  (2.15),  the  Early  voltage  Va  can  be 

expressed  as 


Xg  = - Wepi  + 


V 


Na  (I  + Jc/^a).  qNA(l+Jc/JA)  (2.15) 


v*  = 3! 


(2.16) 
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2.2.3.  Case  3 

Further  increments  in  Jc  drive  the  BJT  into  non-ohmic 
quasi-saturation  (NOQS) . Figure  2. 2d  shows  a typical  electron 
concentration  and  electric  field  distribution  in  NOQS.  In 
NOQS  operation,  the  epi-SCR  excess  f ree-electron 
concentration  (nc  - Nepi)  is  balanced  by  the  heavily  doped 
buried-layer  charge.  Solution  of  Poisson's  equation  leads  to 


XEPI 


26s(Vcb  + foe) 
QNepxCJc/  Jepi  ~l) 


(2.17) 


where  xepi  is  the  epi-layer  SCR  width  in  NOQS  operation  (Fig. 
2. 2d).  In  this  case  Wb  must  be  redefined  as 


Wb  = Wjc  + Wepi  - xepi  • 


(2.18) 


Combining  this  expression  with  (2.10)  and  (2.18)  yields  an 
Early  voltage  expression  for  NOQS  mode: 


w 


b(vCb  + foe) 


0 . 5 Xepi 


KDn  , Esvs(vcb  + foc)jc 

1 4- H 

vs^b  Xepi Wb(Jc  - Jepi)2. 


(2.19) 


The  upper  limit  of  validity  for  (2.16)  and  the  lower  limit 
for  (2.19)  can  be  found  by  replacing  xepi  with  Wepi  in  (2.17): 

Jc  I xEPI  = wEPx  - Jepi  [1+  (xCo/Wepi)2]  (2.20) 

Figure  2.4  shows  a plot  of  the  Early  voltage  using 
(2.11),  (2.16),  (2.19),  and  using  the  "conventional 


it 
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analysis,  (2.12),  assuming  K = 1 and  Dn  = 17  cm2/V*s.  In 

(2.11) ,  (2.16),  and  (2.19),  K is  usually  greater  than  1, 

especially  for  high  currents,  but  this  has  only  a small 
effect  on  the  resulting  Early  voltage. 

In  Fig.  2.4,  discrepancies  between  the  conventional  and 
improved  models  at  low  currents  are  as  Roulston  predicted. 
The  dependence  of  the  SCR  boundary  position  on  electron 
concentration  has  little  effect  for  low  currents.  However, 
for  Jc  greater  than  about  0.1  Jepi,  this  effect  becomes 
significant  and  Va  increases  rapidly  as  Jc  approaches  Jepi 
until  the  entire  epi-layer  becomes  space-charged  (Case  1) . As 
mentioned  earlier  the  potential  drop  across  the  QN  epi  is 
neglected  in  (2.11).  At  low  currents  this  drop  is  negligible. 
As  the  current  approaches  Jepi»  the  IR  drop  becomes  more 
significant,  increasing  Va  slightly  more  than  predicted  by 

(2.11) .  However,  the  resistance  of  the  QN  epi  itself  drops 
with  increasing  current  due  to  the  reduced  width  of  QN  epi. 
Thus  Va  should  again  be  as  predicted  by  (2.11)  at  the 

boundary  between  Case  1 and  2 , where  the  QN  epi  width  reaches 
zero . 

For  Jc  greater  than  Jc  I x = w (Case  2),  the  Early 

voltage  continues  to  increase,  but  at  a slower  rate  than  in 
Case  1.  In  this  case,  the  heavier  buried-layer  doping  reduces 
the  rate  of  base  push-out  with  increasing  Jc- 


EARLY  VOLTAGE  VA  (V) 
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Fig.  2.4.  Comparison  of  conventional  and  improved  VA  models 

for  a transistor  with  NA=lxl017  cm-3,  Nepi=1x1015  cm~ 
3,  Wjc=0.5  |im,  Wepi=7  p.m,  and  Vcb=10  V as  a function 
of  Jc. 
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As  can  be  seen  in  the  figure,  Va  drops  at  the  transition 

between  Case  2 and  3 as  the  device  enters  non-ohmic  quasi- 
saturation. As  shown,  Va  can  even  drop  below  the  prediction 
of  the  conventional  model.  For  Jc  » Jepi/  Va  becomes 
proportional  to  (Jc)1^2. 

For  very  high  base-collector  reverse-biases  or  very 
narrow  epi  thicknesses  the  epi-layer  can  become  completely 
depleted  even  at  low  currents.  For  this  case,  (2.11)  is  not 
valid,  and  (2.16)  is  valid  for  any  current  up  to  the  value  of 
Jc  given  by  (2.20).  Equation  (2.19)  remains  valid  for  NOQS 

operation.  Figure  2.5  shows  the  Early  voltage  variation  as  a 
function  of  Jc  for  K = 1,  Dn  = 17  cm2/V-s,  and  Vcb  = 3 V for  a 

BJT  with  a 0.35  |im  epi  thickness,  which  is  thin  enough  to 
allow  epi-layer  depletion  at  low  currents.  As  shown  in  the 
figure,  carrier-velocity  saturation  can  still  enhance  the 
output  resistance  considerably  in  forward-active  mode  for  Jc 

> Jepi* 


2.3.  Numerical  Simulations 

A series  of  numerical  simulations  were  performed  using 
PISCES-II  and  MMSPICE  to  verify  that  these  effects  can  indeed 
occur.  In  PISCES-II  the  effect  of  carrier-velocity  saturation 
on  Early  voltage  was  studied  by  performing  simulations  with 
and  without  the  field-dependent  mobility  which  models 
carrier-velocity  saturation.  Several  BJT  structures  were 
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COLLECTOR  CURRENT  DENSITY  Jc  (A/cm2) 


Fig.  2.5.  Plot  of  the  improved  Va  model  for  a transistor  with 

Na-0.5x1018  cm-3,  Nepi  = 1x1016  cm-3,  Wjc=0.1  jim, 
Wepi=0.35  |im,  and  Vcb=3  V.  The  epi-layer  is  depleted 
completely  even  at  low  current. 
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simulated,  and  the  results  were  consistent  with  analysis. 
Figure  2 . 6 shows  a comparison  of  Va  behavior  with  and  without 

field-dependent  mobility  for  a transistor  with  a 1.9  Jim  thick 
epi-layer.  As  can  be  seen  in  Fig.  2.6,  using  the  field- 
dependent  mobility  model,  Va  does  increase  considerably  at 

high  currents,  whereas  with  the  field-dependent  mobility 
turned  off,  Va  increases  only  slightly  (due  to  emitter-side 

base-width  modulation) . In  the  lower  curve  (with  field- 
dependent  mobility  turned  off)  the  Early  voltage  drops  at 
high  currents.  This  effect  can  be  explained  by  the  onset  of 
ohmic  quasi-saturation. 

The  physics-based  BJT  model  in  the  circuit  simulator 
MMSPICE  [Je89b]  includes  a charge-based  model  for  advanced 
high-speed  transistors  based  on  one-dimensional  analysis  of 
carrier  transport.  The  analysis  accounts  for  high-current 
effects  and  impact  ionization.  The  model  covers  all  of  the 
cases  considered  in  the  above  analyses.  Appropriate  charge 
equations  are  applied  to  the  various  regions  of  the  device, 
allowing  for  moving  regional  boundaries.  The  resulting  set  of 
coupled  equations  require  iterative  numerical  solution.  The 
highly  physical  basis  of  this  approach  allows  the  software  to 
be  used  for  device  simulation  as  well  as  for  circuit 
simulation. 

We  have  modified  MMSPICE  to  allow  output  of  internal 
physical  parameters  (e.g.,  xb,  xc,  QN  base-width,  and  QN  epi- 
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Fig.  2.6.  PISCES-II  simulation  of  Early  voltage  variation  as 

a function  of  collector  current.  The  device  was 
similar  to  MOSAIC  III  BJT  [8]  except  the  epi  was 

extended  to  1.9  (im.  Nepi  = 1x1016  cm-3,  Wjc=0.15  |im, 

and  Vpr=10  V. 
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width),  in  addition  to  the  usual  voltage  and  current  outputs. 
Figure  2 . 7 shows  Va  and  xb  variations  for  a typical  analog  BJT 
as  a function  of  Ic  for  a fixed  Vce*  As  can  be  seen  in  the 
figure,  MMSPICE  predicts  substantial  increases  in  Early 
voltage  due  to  carrier  velocity  saturation  at  high  currents. 
Figure  2.7  also  clearly  shows  that  the  Early  voltage 
enhancements  are  associated  with  variations  in  xb- 


2.4,  Conclusions 

The  analysis  in  Section  2.2  demonstrates  that  BJT  output 
resistance  can  be  raised  significantly  by  carrier-velocity 
saturation.  For  typical  base-widths,  there  are  usually  small 
increases  in  output  resistance  at  low  currents,  but  the  most 
important  effects  are  predicted  for  currents  close  to  the 
critical  value  Jepi-  These  effects  are  also  predicted  by 
PISCES  and  by  the  highly  physical  circuit/device  simulator 
MMSPICE. 

However,  the  predicted  peak  in  Va  cannot  be  observed  in 
low-frequency  measurements.  The  reason  for  this  inaccuracy  is 
the  neglect  of  the  rise  in  emitter  junction  temperature 
caused  by  an  increase  in  collector  junction  power 
dissipation.  Thermal  feedback  tends  to  severely  degrade  the 
output  resistance  at  high  currents  [Mu64,  Me77,  Fo91a,  Fo91b, 
Fo91c]  which  tends  to  mask  the  Va  enhancement.  With  an 
operating  point  chosen  near  Jepi  it  may  be  possible  to  design 
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high-speed,  high-power  amplifiers  with  very  high  voltage  gain 
if  operation  is  restricted  to  frequencies  high  enough  to 
avoid  thermal  feedback. 
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Fig.  2.7.  MMSPICE  simulation  results  of  Early  voltage  and  SCR 

penetration  into  the  base  as  a function  of 
collector  current  for  a typical  analog  BJT  with 

Nepi  = lxlO15  cm-3,  base-width  = 0.8  |lm,  epi-width  = 

7 Jim,  and  Vcb  = 15  V. 


CHAPTER  3 

SCALABLE  SMALL-SIGNAL  MODEL  FOR  BJT  SELF-HEATING 

3.1.  Introduction 

It  has  long  been  known  [Sp58,  Me77]  that  collector-to- 
emitter  thermal  feedback  can  profoundly  affect  bipolar 
transistor  performance.  This  effect  is  becoming  more 
important  as  current  densities  and  thermal  spreading 
impedances  rise  with  technology  scaling.  Previously  developed 
thermal  impedance  models  were  based  in  the  time  domain 
[ Jo70 ] . It  was  shown  in  [De89]  that  neglecting  thermal 
effects  caused  errors  up  to  4%  in  the  predicted  delay  of  a 
simple  BJT  logic  gate.  This  chapter  will  show  that  local 
heating  causes  much  larger  errors  in  BJT  small-signal 
parameters.  A new  model  is  derived  that  predicts  thermal 
corrections  to  BJT  small-signal  parameters  based  on  a 
solution  of  the  heat-flow  equation  in  the  frequency  domain. 
The  model  is  applicable  in  principle  to  any  emitter  geometry. 
Polynomial  approximations  provide  simplified  expressions  for 
rectangular  emitters.  These  expressions  are  simple  enough  for 
compact  modeling  in  SPICE. 
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3.2.  Model  for  Thermal  Spreading  Resistance 
Figure  3.1  shows  typical  measurements  of  the  emitter 
temperature  response  to  a negative  step  in  collector  power. 
The  measurement  set-up  is  shown  in  Fig.  3.2.  In  Fig.  3.2, 
heating  current  Ih  is  turned  off  during  brief  intervals, 
during  which  Vbe ( t ) is  measured  for  small  fixed  current  Im- 
Electrical  transients  dominate  the  earliest  time  scale  of  the 
measurement,  so  the  temperature  is  plotted  only  for  t > 0.4 
|ls . This  presumably  hides  a short  delay  before  the  emitter 
temperature  begins  to  respond.  The  emitter  temperature 
continues  to  fall  for  several  microseconds.  After  100  |is  the 
local-heating  transient  is  complete  and  the  temperature  is 
constant.  Then,  about  10  ms  after  the  power  step,  the  whole 
package  begins  to  cool;  this  transient  continues  for  many 
minutes.  The  cooling  of  the  package  is  controlled  by  the 
package  thermal  impedance,  which  causes  a coupling  of  all  the 
devices  on  the  chip.  Figure  3.3  shows  the  effect  of  a 
negative  power  step  in  one  transistor's  power  on  the  emitter 
temperature  of  a nearby  transistor.  The  measurement  setup  is 
shown  in  Fig.  3.4.  Only  the  longer  time-scale  transient, 
controlled  by  the  package  thermal  impedance,  is  visible. 
Consistent  with  the  model  of  [Jo70],  thermal  spreading 
impedance  falls  off  strongly  with  distance,  so  that  for 
modest  powers,  even  adjacent  devices  are  thermally  decoupled 


Normalized  Temperature  Change  (°C/W) 
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Time  (s) 


Measured  response  to  negative  collector-power 
step.  Results  of  several  measurements  are 
superimposed. 


Fig.  3.1. 
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Fig.  3.2.  Experimental  setup  for  self -heating . Heating 

current  Ih  is  turned  off  during  brief  intervals, 
during  which  VeE^t)  is  measured  for  small  fixed 
current  Im- 


Normalized  Temperature  Change  (°C/W) 
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Times  (s) 


Fig . 3.3. 


Measured  response  of  transistor  heated  by  nearby 
transistor. 
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i 


Fig.  3.4.  Setup  for  measuring  thermal  response  of  transistor 

heated  by  nearby  transistor. 
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over  the  shorter  time  scale.  This  work  is  concerned  primarily 
with  local  heating  effects. 

If  a unit  impulse  of  heat  occurs  instantaneously  at  some 
point  in  a uniform  medium,  the  temperature  rise  at  a point  a 
distance  r away  is  given  by  [Ca62] 

T (r , t ) = — ~ exp(-^— ) (3.1) 

8pc(7tKt)3/2  4Kt 

where  p is  the  density,  c is  the  specific  heat,  and  K is 

thermal  diffusivity.  An  identical  image  source  can  be  added 
to  create  an  adiabatic  surface  at  the  plane  equidistant 
between  the  sources.  In  their  time-domain  derivation,  Joy  and 
Schlig  [Jo70]  integrated  (3.1)  over  the  volume  of  the 
collector  space-charge  region  (SCR)  and  its  image.  The  result 
is  a closed- form  expression  for  the  impulse  response  at  any 
point  in  the  semiconductor.  The  time  integral  of  this 
expression  gives  the  response  to  a unit  step  in  the  collector 
power.  The  integral  is  complicated  and  does  not  have  a 
closed-form  solution.  The  step  response  at  t = « gives  the 

thermal  resistance  Rth/  defined,  as  the  total  temperature  rise 
caused  by  a unit  step  in  the  collector  power. 

Note  that  (3.1)  is  equivalent  to  thermal  resistance 
which  is  defined  as  the  increments  in  temperature  for  a given 
incremental  unit  power.  Therefore,  a frequency-domain  model 
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can  be  derived  by  first  taking  the  Laplace  transform  of 
(3.1) , yielding 


where  K = pKc  is  thermal  conductivity.  For  s = 0,  this 

equation  gives  Rth  for  the  point  source.  Rth  can  be  integrated 
over  the  collector  SCR  and  its  image  to  yield  the  thermal 
resistance  at  any  point  r'  = (x‘ , y‘ , z')  for  arbitrary 

geometry.  In  rectangular  coordinates 


and  length  of  the  emitter,  D is  the  depth  of  the  base- 
collector  junction,  and  H is  the  SCR  width.  This  integral  is 
not  solvable  in  closed  form,  so  least-squares  fitting  was 
used  to  generate  the  following  approximation  to  (3.3),  in 
which  r'  is  taken  as  the  point  at  the  surface  above  a corner 
of  the  emitter: 

Rth  = l/27CKreff  (3.4a) 

ref f = 2^WL  f if 2 (3.4b) 


47tKr 


(3.2) 


(3.3) 


where  r = V (x-x1 ) 2 + (y-y ' ) 2 + ( z - z ' ) 2 , W and  L are  the  width 


f l (d,h) 


(0.058d+0.14)h+0.34d+0.28 


(3.4c) 
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f2  (a)  = 0.98+0.043a-(6.9-10-4)a1 2+(3.9-10-6)a3  (3.4d) 

where  d = D/Vwl,  h = H/Vwl,  and  a = W/L.  In  (3.4a),  the  2 
instead  of  4 in  the  numetrator  is  given  by  the  fact  that  both 
collector- junction  SCR  and  image  are  considered  as  a point- 
source  of  its  own.  The  approximation  agrees  with  numerical 
solution  of  (3.3)  to  within  5%  for  all  reasonable  geometries. 
Figures  3.5  and  3.6  show  how  fi  and  f2  vary  with  device 
geometry.  Figure  3.7  is  a contour  plot  of  predicted  Rth 

variation  with  position  within  the  emitter  for  a transistor 
with  W = 64  |lm,  L = 48  |i.m,  H = 4 |lm,  and  D = 3 |J.m.  Only  one 

quadrant  is  shown  and  the  center  of  the  emitter  is  at  the 
upper  left.  To  find  Rth  for  any  other  point  r'  at  the  surface 

above  the  emitter,  the  transistor  can  be  broken  into  four 
smaller  rectangular  transistors,  meeting  at  r'.  Equation 
(3.4)  can  be  used  with  modified  values  for  W and  L to  compute 
Rth  for  each  quadrant.  The  overall  value  for  Rth  can  be  found 

using 

1 £ 

Rth  = ^ RthiWiLi  (3.5) 

i=l 

Table  1 gives  calculated  Rth  values  for  r'  located  in  the 
center  of  the  emitter  and  at  the  corner  for  several 
transistor  geometries.  Values  for  H were  computed  using  the 
step- junction  approximation.  We  believe  this  is  the  first 
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Fig.  3.5.  Variation  of  function  fi(d,h)  with  d = D/Vwl,  with 

h = h/Vwl  as  a parameter. 
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Fig.  3.6.  Variation  of  function  f2(a)  with  aspect  ratio 

W/L. 
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Fig.  3.7.  Contour  plot  of  thermal  resistance  at  the  surface 

versus  position  within  emitter  for  a transistor 
with  W = 32  Jim,  L = 24  |lm,  H = 4 |lm,  D = 3 |lm. 
Only  on  quadrant  is  shown  and  the  center  of  the 
emitter  is  at  the  upper  left. 
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compact  model  which  can  predict  the  variation  of  Rth  across 
the  emitter  and  the  variation  with  D and  H. 


3.3.  Small-Signal  Model 

Mueller  [Mu64]  showed  that  BJT  common-emitter  y- 
parameters  can  be  corrected  for  thermal  feedback  using 


YjanE  + Dm^thlmln 


1 - DmZthP 


(3.6) 


where  m or  n = 1 for  the  base  or  2 for  the  collector,  YmnE  is 
the  uncorrected  electrical  Y parameter,  Zth  is  the  thermal 
impedance,  and  P is  the  dissipated  power.  D2  is  the 

fractional  temperature  coefficient  of  the  collector  current, 
(typically  about  7 %/K)  and  Di  is  that  of  the  base  current, 
equal  to  D2  - Dp,  where  Dp  is  the  fractional  temperature 
coefficient  of  p,  (typically  around  0.7  %/K).  Equation  (3.6) 
can  be  used  with  Zth  = Rth  to  thermally  correct  the  DC  small- 
signal  parameter  gmn. 

Figure  3.8  compares  measured  and  simulated  values  for 
the  normalized  output  resistance  parameter  Va  = Ic/<?22  ~ Vce 
for  a silicon  BJT  with  23  x 23  (|J.m)2  emitter,  collector 
junction  depth  of  3.0  fim,  and  epi-doping  of  1015  /cm3.  The 
value  for  H in  (3.4)  was  2.8  |im,  calculated  assuming  a step 
junction  with  Vce  = 6 V,  which  led  to  Rth  = 70  °C/W.  The 
simulation  used  g22E  = Ic/Vae#  where  Vae  = 350  V is  the 


Output  Resistance  Parameter 
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Collector  Current  (A) 


Fig . 3.8. 
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Mill  Hi 

1CT2 


Measured  and  simulated  normalized  output 
resistance  Va  = lc/922  - Vce-  Simulations  are  based 
on  (3.4)  and  (3.6) . 
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Table  1:  Calculated  and  measured  thermal  resistance  for 

several  transistor  geometries.  Calculations  were 
based  on  (3.4)  and  (3.5).  Measured  values  give  best 
fits  to  output  resistance  and  reverse 
transconductance  using  (3.6)  and  include  an 
estimate  of  experimental  error. 


Type 

W (|!m) 

L (|Xm) 

Rth  at 
corner 

Rth  at 
center 

Measured 

Rth 

23 

23 

69 

115 

74+5 

a 

83 

23 

36 

65 

41±4 

12 

12 

112 

170 

118±10 

4 

2 

492 

778 

480+20 

b 

8 

2 

348 

582 

3 90±5 

10 

2 

306 

520 

350±20 

10 

4 

251 

443 

290±10 

c 

7.2 

1.2 

467 

814 

435±20 

a: 

D = 

3 |lm,  Nepi  = 

lxlO15  cm 

~3  / Vce  = 6 

V,  H = 2.83  |J.m. 

b: 

D = 

0.4  |lm,  Nepi 

= 7xl015 

cm-3 , Vce  = 

3 V,  H = 0 

.85  |0,m. 

c: 

D = 

0.25  [lm,  Nepi 

= 2xl016 

cm-3,  Vce 

= 2 V,  H = 

0.35  |lm . 
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electrical  Early  voltage,  measured  at  low  currents.  The  data 
were  corrected  to  account  for  package  thermal  resistance  of 
110  °C/W,  measured  using  the  plateau  near  1 ms  shown  in  Fig. 
3.1.  Note  that  most  circuit  models  use  a constant  value  (the 
Early  voltage)  for  this  parameter.  Measured  g22  or  912 
variations  can  be  used  in  (3.6)  to  extract  a best-fit  value 
for  Rth*  The  last  column  in  Table  1 compares  Rth  values  found 
by  these  methods  to  calculated  values.  This  method  will  be 
discussed  in  Chapter  4.  Note  that  the  measured  values  match 
well  with  the  calculated  Rth  for  the  corner  of  the  emitter. 

This  result  appears  to  be  largely  coincidental,  but  it  is 
fortunate,  as  it  means  that  a single  effective  temperature 
can  be  assigned  to  the  emitter  using  (3.4)  with  good 
accuracy. 

The  other  DC  small-signal  parameter  which  is  strongly 
affected  by  thermal  feedback  is  gi2-  Electrical  models  predict 

only  a small  negative  value  for  this  parameter,  caused  by 
Early-effect  modulation  of  recombination  in  the  quasi-neutral 
base.  This  effect  is  normally  neglected.  Equation  (3.6) 
predicts  a substantial  positive  value  for  gi2  approaching 
g22/P  with  increasing  Ic.  This  can  have  a substantial  effect 

on  circuit  performance,  as  the  input  and  output  impedances 
become  strongly  dependent  on  the  loading  at  the  opposite 
port.  Note  that  (3.6)  shows  that  errors  in  small-signal 
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parameters  due  to  neglecting  thermal  feedback  can  occur 
without  large  DC  power  dissipation. 


3.4.  AC  Modeling 

Equation  (3.2)  can  be  used  to  model  thermal  spreading 
impedance  above  DC  if  the  reff  from  (3.4)  is  substituted  for  r 
in  (3.2).  The  physical  significance  of  reff  is  an  effective 

distance  between  the  collector- junction  SCR  and  the  point  at 
the  surface  above  a corner  of  the  emitter.  The  results  are 
compared  to  the  single-pole  model  [De89,  Vo89,  Ye91]  and  the 
inverse  Fourier  transform  of  Joy  & Schlig's  predicted  impulse 
response  in  Fig.  3.9.  For  this  case  (the  same  device  as  in 
Fig.  3.1)  the  magnitude  of  the  point-source  model  differs  by 
less  than  6%  compared  to  the  Fourier  transform  of  the  full 
time-domain  model  for  frequencies  up  to  fpk  =350  kHz,  defined 
as  l/(27ttpk),  where  tpk  = reff2/(6K)  = 0.45  (is  is  the  time  at 

which  the  impulse  response  in  (3.1)  reaches  its  maximum.  At 
this  frequency  iZthl  has  dropped  to  18%  of  Rth«  Note  that  the 
commonly  used  one-pole  approximation  predicts  most  of  the 
drop  in  IZthl  to  occur  over  one  decade  of  frequency,  whereas 
the  IZthl  in  Fig.  3.9  varies  over  three  decades  of  frequency. 

The  phase  response  is  less  accurate,  differing  by  more 
than  30  degrees  for  frequencies  greater  than  fpk.  These  phase 

errors  are  easily  understood:  the  model  represents  the 
transistor  as  a point-source  collector  and  a point  emitter, 


Phase  of  Zth  ( ) Magnitude  of  Zth  (°C/W) 
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Frequency  (Hz) 
(a) 


F r equency  (Hz) 
(b) 


Fig.  3.9.  Magnitude  and  phase  of  thermal  spreading  impedance 

Zth-  Solid  lines  were  computed  using  (3.2)  and 
(3.4) . 
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separated  by  an  reff  which  is  inevitably  greater  than  the  base 

width.  Thus,  the  model  overestimates  the  initial  delay. 
However,  since  they  coincide  with  decreasing  I Zth  I » these 
phase  errors  cause  relatively  little  error  in  the  y- 
parameters.  Figure  3.10  shows  typical  measured  and  simulated 
values  for  Y22 / both  magnitude  and  phase  for  Vbe  = 0.8  V,  Vce 
= 6 V,  and  Iq  = 3.74  mA.  As  can  be  seen  in  the  figure, 

simulations  show  good  agreement  with  measurements  despite  the 
inaccuracy  in  the  phase  of  the  model.  Figures  3.9  and  3.10 
prove  the  importance  of  correctly  doing  thermal  modeling  for 
BJT  output  conductance. 

3.5.  Conclusion 

The  model  described  above  provides  a simple  and 
practical  way  to  improve  the  accuracy  of  BJT  circuit 
simulation.  Polynomial  approximations  are  used  to  express  the 
variation  of  Rth  across  the  emitter;  the  results  are 
equivalent  to  the  infinite  integral  in  [Jo70] . Measurements 
show  that  a single  value  of  Rth  can  represent  typical  BJT 
thermal  behavior  at  DC.  A frequency-domain  solution  for  Zth  is 
presented  which  is  more  accurate  than  the  commonly  used 
single-pole  thermal -impedance  model. 

Considering  the  significant  effects  which  thermal 
feedback  has  on  BJT  behavior,  it  is  interesting  to  speculate 
why  thermal  effects  are  not  yet  routinely  included  in  circuit 
simulators,  in  spite  of  prior  work  on  including  such  effects. 
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Frequency  (Hz) 


(a) 


Frequency  (Hz) 

(b) 

Fig.  3.10.  Measured  and  simulated  Y22  for  Vbe  = 0-8  V,  Vce  = 6 

V,  and  Ic  = 3.74  mA.  (a)  Magnitude,  (b)  Phase. 
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[Fu76,  Vo89] . One  contributing  factor  may  be  the  willingness 
of  many  circuit  designers  to  accept  empirical  models  --  model 
parameters  can  be  adjusted  to  account  for  many  thermal 
effects.  Figure  3.11  demonstrates  that  even  this  empirical 
approach  can  lead  to  significant  errors.  Extractions  of  SPICE 
parameter  VAF  were  done  under  several  different  conditions. 
Differences  greater  than  30%  were  found  depending  on  HP4145 
integration  time,  and  whether  or  not  a fan  was  blowing  on  the 
package  to  change  its  effective  thermal  impedance.  This 
should  give  pause  to  the  analog  IC  designer  --  such  thermal 
effects  could  cause  the  gain  of  a two-stage  op  amp  to  vary 
more  than  2 to  1 depending  on  package  type  and  mounting 
conditions . 


Early  Voltage 
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Collector  Current  Ic  [A] 


Fig.  3.11.  Vaf  values  optimized  for  various  currents, 

extracted  under  different  thermal  conditions:  + 

using  an  HP-4145  at  minimum  delay  between  points 
*:  static  measurement  with  small  fan  blowing  on 

the  14-pin  DIP;  o:  static  measurement  without  fan. 
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CHAPTER  4 

PARAMETER  EXTRACTION  FOR  BIPOLAR  CIRCUIT  MODELING 

4.1  Introduction 

Thermal  feedback  from  collector  power  to  emitter 
temperature  has  been  shown  to  significantly  affect  BJT 
performance  [Mu64,  Jo70,  Me77,  De89].  Recently  several 

workers  [Vo89,  Ye91]  have  demonstrated  ways  to  model  thermal 
feedback  in  circuit  simulation  assuming  thermal  model 
parameters  are  available.  This  chapter  introduces  simple  ways 
to  extract  thermal  parameters  needed  for  such  simulations. 

BJT  thermal  behavior  is  dominated  by  two  components, 
characterized  by  the  thermal  spreading  impedance  and  the 
chip/package  thermal  impedance.  It  is  important  to  extract 
these  components  separately,  as  the  package  thermal  impedance 
causes  a coupling  between  transistors  on  the  same  chip,  and 
thermal  spreading  impedance  usually  does  not.  Standard 
methods  for  extracting  thermal  impedance  are  adapted  for 
characterizing  power  transistors  [Oe76].  They  are  poorly 
suited  for  measuring  thermal  spreading  impedance,  which 
becomes  increasingly  important  with  decreasing  transistor 
size . 
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4.2.  Thermal  Impedance  Extraction  Methods 
The  standard  electrical  method  for  thermal  impedance 
measurement  is  the  "emitter-only  switching"  method.  In  this 
technique,  a large  heating  current  is  first  applied  to 
establish  thermal  equilibrium.  This  current  is  then  turned 
off  and  the  temperature  is  inferred  from  the  resulting 
transient  response  of  Vbe/  measured  at  a small  fixed  emitter 
current.  The  initial  response  is  controlled  by  both 
electrical  and  thermal  spreading-impedance  effects.  This 
makes  it  difficult  to  determine  the  original  temperature,  and 
prevents  accurate  extraction  of  the  thermal  spreading 
resistance  Rth-  The  method  is  still  useful,  however,  for 
finding  package  thermal  impedance.  For  small  transistors,  the 
transient  due  to  thermal  spreading  impedance  is  typically 
completed  long  before  the  thermal  transient  of  the  package 
begins.  The  package's  transient  typically  begins  a few 
milliseconds  after  the  power  step  and  lasts  for  up  to  several 
minutes  [Fo91a] . The  package  thermal  resistance  can  be  found 
by  dividing  the  total  change  in  temperature  after  1 ms  by  the 
size  of  the  power  step. 
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4.2.1.  Using  output  resistance  as  a temperature sensitive 

parameter 

The  following  technique  can  be  used  to  extract  thermal 
spreading  resistance.  It  is  easy  to  show  [Mu64]  that  for 
small  changes  in  the  base  and  collector  voltages,  the  change 
in  Ic  can  be  written  as 


Ale  = <?21eAVbe  + <?22eAVce  + IcDCRth  ( IcAVce  + Vce^Ic)  (4.1) 

where  gi2E  and  g22E  are  common- emit ter  small-signal  two-port 
parameters  neglecting  thermal  effects,  and  Dc  is  the 
fractional  temperature  coefficient  of  collector  current  (Dc  = 
Ic_1  3lc/3T)  . With  AVBe  set  to  zero,  (4.1)  can  be  solved  for 
DcRth»  which  gives 


s 


DCRth 


Alc/Ic  ~ AVce/  (Va  + Vce) 
Ic^Vce  + Vce^Ic 


(4.2) 


where  g22E  is  expressed  as  Tc/  (Va  + Vce)  » where  Va  is  the 
Early  voltage. 

The  extraction  procedure  is  as  follows:  Vce  is  set  to 

some  value  in  the  forward-active  region.  For  a series  of 
fixed  values  of  Vbe»  Vce  is  varied  slightly  and  the  resulting 
Ale  is  measured.  Since  thermal  effects  are  negligible  for  low 
currents,  a constant  value  for  Va  can  be  extracted  at  low  Vbe 
values  using 


VA  = (AVce/AIc)  ic  - VCE 


(4.3) 
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At  higher  currents,  thermal  effects  become  significant,  and 
(4.2)  can  be  used  to  extract  the  product  DcRth-  For  moderate 
currents,  Dc  can  be  approximated  as  [Me77] 


1 

T 


Vgq  - Vbe 
L VT 


(4.4) 


where  T is  absolute  temperature,  Vt  is  the  thermal  voltage, 
Vqo  for  silicon  is  1.206  V,  and  T|  = 2 (related  to  the  doping 

level).  Equation  (4.4)  becomes  inaccurate  at  very  high 
currents,  so  it  is  important  not  to  use  such  high  currents  in 
the  extraction. 

It  is  important  to  carefully  consider  the  timing  of 
these  measurements  to  avoid  confusing  thermal  spreading  and 
package  thermal  impedances.  An  HP-4145  Semiconductor 
Parameter  Analyzer  set  to  its  minimum  integration  time 
requires  at  least  20  ms  to  measure  a pair  of  (Vce <■  Ic)  data 
points,  during  which  time  substantial  change  in  package 
temperature  can  occur.  Longer  delays  between  data  points  lead 
to  larger  errors . (For  this  reason  faster  equipment  is  needed 
to  extract  thermal  spreading  resistance  alone?)  A pulse 
generator,  a current-to-voltage  converter  and  an  oscilloscope 
are  generally  adequate.  Figure  4.1  show  the  measurement  set- 
up using  such  equipment.  In  the  figure,  the  pulse  generator 
varies  the  collector  current  of  the  device  for  a given  Vbe- 

Variations  in  collector  current  are  amplified  by  the  op-amp 
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Fig.  4.1  Pulse  measurement.  (a)  Experimental  setup  to 

measure  Rth  using  a pulse  generator,  a current-to- 
voltage  converter,  and  an  oscilloscope,  (b)  A 
typical  shape  of  measured  output. 
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and  measured  with  an  oscilloscope.  Variation  in  collector 
current  is  given  by  Ale  = Av/R,  where  Av  is  the  peak-to-peak 

voltage  of  op-amp  output  voltage  and  R is  the  op  amp's 
feedback  resistor.  In  Fig.  4.1,  the  pulse  width  tw  has  to  be 

short  enough  to  avoid  any  error  caused  by  package  thermal 
impedance . 

Another  method  which  avoids  the  effect  of  package 
thermal  impedance  involves  measuring  the  output  admittance  Y22 

using  an  impedance  analyzer.  As  previously  noted,  for  typical 
small  transistors  there  is  a period  after  a power  step  during 
which  the  thermal  step  response  is  constant  for  some  time. 
This  implies  that  there  is  also  a range  of  frequencies  over 
which  the  thermal  frequency  response  is  constant.  In  the 
middle  of  this  frequency  range  Y22  is  also  constant  for  a 
given  current.  Its  real  part  can  be  denoted  g22  • At  low 
currents,  thermal  effects  are  negligible  and  g22  = g 22E  = 
lc/  (Va  + Vqe)  / which  can  be  used  to  find  Va-  With  increasing 
current,  g22  rises  above  Iq/  (Va  + Vce)  > and  Rth  can  be 
extracted  using 

_ 1 g22/ic  - 1/ (Va  + Vce)  ia  ^ 

Rth  = 5^  lc  - g22vCE  <4'5) 

This  method  generally  gives  results  consistent  with  those 
found  with  the  previous  method.  Figure  4.2  shows  a typical 
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plot  of  Rth  extracted  at  various  collector  currents  based  on 
the  above  methods. 

4.2.2.  Using  reverse  transconductance  as  a temperature 

sensitive  parameter 

In  the  characterization  method  described  above,  the 
Early  voltage,  Va  has  to  be  supplied  to  extract  Rth  from  a 
measured  g22  • For  typical  analog  devices,  the  effect  of 
emitter-side  base-width  modulation  ("late  effect")  is  usually 
negligible  so  that  the  Early  voltage  stays  roughly  constant 
with  variation  in  base-emitter  voltage  for  a given  collector- 
base  voltage.  In  this  case,  the  Early  voltage  measured  at  low 
currents  can  be  used,  with  little  error,  to  extract  Rth  from 
the  g22  measured  at  other  currents.  However,  for  typical 
advanced  BJTs,  late  effects  are  far  more  significant  because 
of  the  narrower  base  width,  so  that  the  extracted  Early 
voltage  varies  considerably  as  a function  of  base-emitter 
voltage.  In  this  case,  using  the  Early  voltage  measured  at 
low  current  can  leads  to  significant  error  in  extracted  Rth* 
Extracting  Early  voltage  at  high  currents  is  not  easy  because 
of  the  thermal  effects.  It  may  still  be  possible  to  extract 
Rth  using  the  previous  method,  but  can  be  cumbersome,  and  not 


Thermal  Resistance  (°C/W) 


61 


Collector  Current  (A) 


Fig.  4.2 


A typical  plot  of  Rth  extracted  at  various  Ic. 
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very  precise.  A better  method  for  this  case  uses  the  reverse 
transconductance  gi2  as  the  temperature-sensitive  parameter. 
Recently,  Reisch  [Re92]  also  have  independently  suggested 
this  method  that  will  be  described  in  the  following. 

An  equation  comparable  to  (4.1)  gives  the  change  in  base 
current  as 

AIb  = CJIIeAVbe  + 912E^Vce  + iB^Rth  ( IcAVce  + VceAIc)  (4.6) 

where  Db  is  the  fractional  temperature  coefficient  of  base- 
current  Ib*  For  the  typical  advanced  devices,  gi2E  is 
negligible.  With  AVbe  set  to  0,  (4.6)  can  be  solved  for  DBRth/ 
which  gives 


With  similar  precautions,  the  equipment  described  in  the 
previous  extraction  method  can  be  applied  for  this  method. 
For  moderate  currents,  assuming  base  currents  are  dominated 
by  back  injection  into  the  emitter  [Je91b]  , Db  can  be 

approximated  as  [S176,  Sz81] 


E>BRth  = 


AlB 


(4.7) 


Ib ( IcAVce  + VceAIc ) 


db  = ~ [ 


1 .Vgq  ~ VbE  - AVg 


+ ri] 


(4.8) 


VT 


where  AVq  is  the  bandgap  narrowing  which  is  a function  of  the 

impurity  doping  concentration  in  the  emitter.  The  accuracy  of 
the  extracted  Rth  depends  on  the  accuracy  of  AVg  and  T| . 
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However,  the  errors  are  not  significant  as  these  parameters 
are  typically  small  compared  to  the  other  terms.  Equation 
(4.8)  becomes  inaccurate  at  very  high  currents  due  to  the 
base  and  emitter  parasitic  resistances. 

Note  that  output  resistance  should  always  be  extracted 
with  the  base  voltage  fixed,  not  the  base  current.  Since  gi2E 
is  extremely  small  for  narrow-base  transistors,  fixing  Ib  and 
fixing  Vbe  would  be  equivalent  in  the  absence  of  thermal 
effects.  In  practice,  however,  with  fixed  T-q,  Vbe  adjusts 
itself  to  cancel  most  of  the  thermal  effect  on  the  output 
resistance.  This  effect  is  easily  demonstrated  by  comparing 
output  characteristics  measured  with  fixed  Ib  and  with  fixed 
Vbe-  At  high  currents,  the  fixed-VBE  curves  generally  show 
much  greater  slopes. 


4.3.  Extraction  Results 

Measured  thermal  spreading  resistance  values  were  used 
with  (4.1)  and  (4.6)  to  predict  BJT  small-signal  behavior. 
Values  for  Dc  and  Db  were  determined  using  SPICE.  Figure  4.3 
shows  measured  and  simulated  variation  of  normalized  output 
resistance  parameter  Va'  = rQ  Ic  - Vce  with  Ic-  The  transistor 
had  an  emitter  area  of  23  x 23  (p.m)2,  base  width  of  1.0  fim, 
and  epitaxial-layer  doping  of  1015  cm-3.  The  measured  thermal 
spreading  resistance  Rth  was  70  °C/W.  The  device  was  packaged 

in  a 14-pin  plastic  DIP  whose  chip-to-ambient  thermal 
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Fig.  4.3  Measured  and  simulated  normalized  output  resistance 

parameter  Va  = Ic/Q'22  - Vce  for  a 23  x 23  (|lm)2  BJT . 
Rth  = 70  °C/W. 
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resistance  was  110  °C/W.  The  measurements  and  the  simulations 
were  both  arranged  to  include  only  the  thermal  spreading 
resistance.  Figure  4.4  shows  measured  and  simulated  variation 
of  gi2  for  the  same  transistor.  In  this  case,  an  HP4145 

Parameter  Analyzer  set  to  its  maximum  integration  time  was 
used  to  make  the  measurement.  As  discussed  previously,  this 
causes  the  package  thermal  impedance  to  contribute  a 
component  to  the  effective  value  of  Rth-  Rth  was  raised  to  90 
°C/W  in  the  simulation  to  account  for  this  effect. 

The  measurements  and  simulations  in  both  figures  agree 
well  over  many  decades  of  current.  The  disagreement  at  high 
currents  could  be  due  to  effects  of  parasitic  emitter  and 
base  resistances  or  errors  in  the  models  for  Dc  and  Db  at 
these  current  levels.  Note  that  SPICE  models  predicts  a 
constant  value  (Vaf)  for  normalized  output  resistance 
parameter  and  a zero  or  very  small  negative  value  for  gi2-  The 
much  larger  positive  value  for  gi2  seen  here  is  important 
because  it  makes  the  input  impedances  depend  strongly  on 
loading  condition  at  the  opposite  port. 

Similar  agreement  between  measured  and  simulated  small- 
signal  parameters  have  been  seen  with  many  other  transistors. 
Table  4.1  shows  Rth  values  found  by  the  characterization 

methods  described  in  this  chapter.  Of  special  interest  is  the 
high  thermal  spreading  resistance  seen  with  dielectrically 
isolated  BJTs.  Such  transistors  are  isolated  from  the  silicon 
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Fig.  4.4  Measured  and  simulated  reverse  transfer  conductance 

gi2  for  23  x 23  (Jim)2  BJT.  Simulation  used  Rth  = 90 
°C/W  to  account  for  package  thermal  impedance 
affecting  the  measurement. 
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Table  4.1  Measured  thermal  spreading  resistance  for  several 

transistor  geometries . 


Number  of 


Type W (|Xm) L (|im) emitter Measured  Rth 


9 

9 

1 

118±10 

9 

9 

2x4 

50 

54 

9 

2 

50 

9 

9 

8 

42 

a 

20 

20 

1 

7 4±5 

20 

20 

2 

57 

80 

20 

1 

50±5 

20 

20 

4 

43 

80 

20 

2 

40±2 

20 

20 

8 

31 

4 

2 

1 

480±20 

b 

8 

2 

1 

390±5 

10 

2 

1 

350120 

10 

4 

1 

290110 

c 

4 

2 

2 

560120 

8 

2 

1 

435120 

d 

10 

6 

1 

280120 

e 

22 

11 

1 

32713 

22 

22 

2 

28015 

f 

22 

22 

2 

350120 

a: 

NPN, 

Base-width 

= 1 |im,  Nepi  = 

lxlO15 

cm-3 . 

b: 

NPN, 

Base-width 

= 0.5  |im,  Nepi 

= 7x10 

15  cm-3 

♦ 

c: 

NPN, 

Base-width 

= 0.15  (Xm,  Nepi  = 2xl016  cm 

3 

♦ 

d: 

NPN, 

similar  vertical  dimension  as  c. 

» 

e : 

NPN, 

DI,  Base-width  = 0.8  |im. 

Nepi  = 

lxlO15 

cm-3 . 

f : 

PNP, 

DI,  Base-width  = 0.8  |im. 

Nepi  = 

lxlO16 

cm-3 . 

- W and  L represent  the  dimensions  of  the  mask  opening. 
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substrate  by  a layer  of  silicon  dioxide,  whose  thermal 
conductivity  is  only  1 % that  of  silicon.  The  Rth  of  a 23  x 23 

(fim)2  dielectrically  isolated  transistor  was  found  to  be  280 
°C/W,  four  times  that  of  the  similar-size  junction-isolated 
transistor  previously  described.  Note  that  the  thermal 
resistance  of  the  vertical  PNP  transistor  (dielectrically 
isolated)  shows  a comparable  but  slightly  higher  value  of  Rth 

compared  to  the  same-dimension  NPN  transistor.  Slightly 
higher  value  of  Rth  for  the  PNP  transistor  is  consistent  with 


the  model  described  in  Chapter  3 (Fig.  3.5)  considering  the 
higher  epi-doping  concentration. 


Modeling  of  thermal  spreading  impedance  is  expected  to 
become  increasingly  important  in  circuit  and  device 
simulation  as  transistor  geometries  are  scaled  down  and 


to  extract  thermal  spreading  resistance  of  BJTs  separately 
from  package  thermal  impedance.  These  techniques  will  be 
especially  useful  for  modeling  of  devices  whose  unusual 
geometries  or  structures  (such  as  dielectric  isolation)  make 
analytical  modeling  of  thermal  impedance  impractical. 


4.4.  Conclusion 


current  densities 


chapter  presents  practical  ways 


CHAPTER  5 

THERMAL  EFFECTS  IN  BJT  INTEGRATED  CIRCUITS 


5.1.  Introduction 

Circuit  and  device  integration  levels  continue  to 

\ 

improve  with  time.  The  resulting  smaller  device  sizes  result 
in  higher  thermal  spreading  resistances  and  the  demand  for 
higher  performance  gives  rise  to  higher  power  densities  in 
these  devices.  A result  of  these  trends  is  towards  higher 
junction  temperature  elevations  in  the  devices.  In  addition, 
in  new  isolation  technologies  (e.g.,  trench  isolation, 
dielectric  isolation,  and  SOI),  the  poor  thermal  conductivity 
of  silicon  dioxide  makes  the  temperature  elevations  far  more 
significant . 

As  discussed  in  the  previous  chapters,  for  modest 
powers,  even  adjacent  devices  are  thermally  decoupled  from 
each  other  [Fo91a]  so  that  localized  heating  becomes 
intrinsic  to  the  BJT  itself.  This  work  focuses  on  this 
intrinsic  self-heating  effects  which  can  not  be  avoided 
through  changing  the  layout  of  the  IC . 

It  was  shown  in  [De89]  that  neglecting  intrinsic  local 
heating  effects  caused  errors  up  to  4 % in  the  predicted 
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delay  of  a simple  BJT  logic  gate.  Local  heating  causes  much 
larger  errors  in  some  widely  used  high  accuracy  monolithic 
linear  integrated  circuits  such  as  D/A  converters, 
instrumentation  amplifiers,  analog  multipliers,  logarithmic 
amplifiers,  current  mirrors,  and  precision  voltage 
references.  Local  heating  can  also  cause  errors  in  op-amp 

v 

practice  as  it  can  lower  the  low-frequency  open-loop  gain, 
cause  distortion,  and  change  the  phase  response.  Therefore, 
it  is  important  for  designers  to  have  good  insight  into  the 
circuit  implications  of  local-heating  effects,  as  even  well- 
executed  designs  can  fall  prey  to  these  effects  [To90] . 

Knowledge  of  the  implications  of  thermal  effects  in 
integrated  circuits  not  only  provides  the  design  principles 
to  reduce  thermal  effects  in  IC  design  but  can  be  used  to 
avoid  them  for  better  IC  design.  The  object  of  this  chapter 

S 

is  to  provide  general  insight  into  local-heating  effects  on 
the  DC  and  small-signal  performance  of  bipolar  integrated 
circuits  and  to  identify  the  circuit  conditions  under  which 
self-heating  causes  significant  errors  in  simulation.  The 
thermal  impedance  model  described  in  Chapter  3 has  been 
implemented  in  SPICE  DC  and  AC  analysis  by  D.  Zweidinger 
[Zw92]  . This  improved  BJT  model  is  used  to  verify  the 
discussions  in  this  chapter. j 

. — . ^ 1 i 

In  Section  5.2,  the  improved  SPICE  BJT  model  [Fo91c]  is 
introduced  which  includes  a brief  description  of  how  the 
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thermal  models  are  incorporated  in  the  DC  and  AC  routines, 
and  discusses  the  three  levels  of  thermal  modeling  available 
to  the  user  [Zw92] . In  addition,  a series  of  simulations  are 
compared  to  measurements  to  demonstrate  the  improved  accuracy 
afforded  by  the  thermally  enhanced  BJT  model. 

In  Section  5.3,  DC  thermal  effects  are  considered  in  two 
aspects,  magnification  of  the  increment  in  collector  current 
due  to  the  mismatch  in  collector-base  voltage  and  collector 
current  deviation  from  near-exact  exponential  dependence  on 
base-emitter  voltage.  For  small-signals,  the  effects  of  self- 
heating on  the  input  and  output  resistances  and  the  effects 
of  series  emitter,  base,  and  collector  resistances  are 
discussed.  The  effects  of  discrepancies  in  the  input  and 
output  impedance  on  other  small-signal  behavior  are  also 
discussed.  Simulation  results  of  a few  circuit 
configurations,  based  on  modified  and  ordinary  SPICE,  are 
provided  for  verification  and  demonstration  purposes.  Design 
strategies  to  avoid  significant  thermal  problems  are  also 
discussed  in  this  section. 

In  Section  5.4,  the  work  discussed  in  this  chapter  are 
summarized. 

5.2.  Modified  SPICE  BJT  model 

From  Chapter  4,  a simple  model  for  the  thermal  impedance 
of  any  rectangular  emitter  geometry  is  given  by 
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Zth  = Rthexp(-reffVs/k) 


where 


reff  = 


1 


2 JtKRth 


and  Rth  is  given  as 


Rth  = 1/4  JIKVWL  fi  f2 


(5.1a) 


(5.1b) 


(5.1c) 


where 

f 1 (d, h)  = (0.058d+0.14)h+0. 34d+0 . 28 , (5. Id) 

f2  (a)  = 0 . 98+0 . 043a- (6 . 9 • 10-4) a2+ (3 . 9 • 10-6)  a3 , (5.1e) 


a = W/L, 


(5. If) 


d = D/VwL, 


(5 . lg) 


h 


2e  (Vp  + VCB) 

QNepi 


( 5 . lh) 


Vc  is  the  collector-base  junction  built-in  potential.  This 
model,  which  is  valid  for  junction-isolated  devices,  has  been 
implemented  in  SPICE  DC  and  AC  analysis.  MicroSim's  PSpice 
[Ps89]  was  used  as  the  initial  code  for  the  modifications, 
using  the  device  models  option.  The  SPICE  DC  analyses  were 
modified  in  three  respects  to  include  thermal  effects.  First, 
in  the  operating-point  routine,  the  large-signal  currents  and 
all  temperature -dependent  parameters  were  adjusted  to  reflect 
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a raised  local  temperature  for  each  device  using  the  standard 
SPICE  equations.  Second,  the  admittances  were  modified  to 
include  thermal  effects  using  (3.6).  Finally,  the  increase  in 
temperature  from  one  iteration  to  the  next  was  limited  to 
improve  convergence  of  the  operating  point  and  prevent 
thermal  runaway.  Figure  5.1  shows  a flowchart  of  the  modified 
subroutine  to  load  the  conductance  matrix,  which  is  the  heart 
of  the  DC  analysis  and  operating-point  computation.  The 
modifications  to  the  AC  analysis  routine  consist  primarily  of 
changes  to  the  admittances,  which  are  similar  to  those  in  the 
DC  analysis. 

The  modified  SPICE  requires  a few  parameters  to  be  added 
to  the  BJT  model  description.  Three  levels  of  thermal 
modeling  are  provided  to  the  user.  Each  level  requires  a 
different  way  of  specifying  thermal  impedance  parameters  in 
the  model  specification.  The  first  level  is  fully  predictive 
and  requires  only  information  about  the  device  geometry  and 
doping.  From  (5.1),  the  thermal  impedance  can  be  specified 
for  a given  W,  L,  D,  and  Nepi.  This  predictive  approach  allows 

circuits  to  be  simulated  concurrently  with  a developing 
process . 

The  second  method  of  determining  Zth  uses  a value  of  Rth 
extracted  from  measurements.  The  frequency  response  is  then 
predicted  using  (5.1a)  and  (5.1b).  Suggested  techniques  to 
measure  Rth  are  given  in  Chapter  4.  This  modeling  approach  can 
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Fig.  5.1.  Flowchart  of  BJT  DC  analysis. 
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give  greater  accuracy  than  the  predictive  approach  when  good 
measurements  of  Rth  can  be  made. 

(The  third  approach  uses  a single-pole  model,  in  which 
Zth  is  specified  by  providing  the  real  part,  Rth,  and  the 
imaginary  part,  Cth-  When  specified  in  the  manner,  simulations 

are  only  accurate  over  a limited  frequency  range. 

A series  of  simulations  of  Ic  vs.  Vce/  Early  voltage  Va 

vs.  Ic,  gi2  vs.  Ic,  and  Y22  vs.  frequency  are  compared  to 
measurements  in  Figs.  5.2  - 5.5,  respectively.  The  transistor 
used  was  rather  big.  Its  thermal  spreading  resistance  Rth  was 
near  70  °C/W.  The  device  had  emitter  length  and  width  of  23 
|im,  collector-base  junction  depth  of  3.0  |im,  and  epi-doping 
Nepi  of  1015  cm-3.  The  fully  predictive  model  was  used  for  the 
simulations.  In  Fig.  5.2,  as  expected,  measurements  and 
modified  SPICE  simulations  predict  higher  collector  current 
and  much  higher  Ic  - Vce  slopes  than  ordinary  SPICE  does.  The 
increase  in  Ic  - Vce  slope  for  high  Ic  can  be  seen  more 
clearly  in  Fig.  5.3.  The  Early  voltage  VA,  which  is  predicted 

to  be  constant  based  on  the  electrical -only  model,  approaches 
zero  for  high  Ic  as  the  thermally  induced  variation  dominates 
the  Ic  - Vce  slope.  Note  that  this  transistor  has  a rather 
high  electrical-only  Early  voltage  Va  so  that  even  though  Rth 

is  only  70  °C/W,  the  output  resistance  drops  to  half  of  what 
the  electrical-only  model  predicts  at  only  Ic  = 400  (1A. 


77 


u 

H 

-U 

a 

a) 

u 

u 

d 

a 

o 

4J 

u 

a) 


o 

u 


Fig.  5.2.  Typical  measured  and  simulated  collector  currents 

as  a function  of  collector-emitter  voltage. 
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Collector  Current  Ic  (A) 


Fig.  5.3.  Normalized  output  resistance  = Ic/922  _ Vce  as  a 

function  of  Ic,  measured  and  simulated.  Vce  was  set 
to  6V.  Note  that  unmodified  Gummel-Poon  model 
predicts  a constant  for  this  parameter. 
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In  Fig.  5.4,  measurement  and  modified  SPICE  simulation 
of  gi2  as  a function  of  Ic  show  good  agreement  up  to 
reasonably  high  currents.  The  discrepancies  at  high  currents 
are  probably  caused  by  the  inaccuracies  in  the  model  for  the 
temperature  coefficient  of  Ib  and  by  inexact  parameters  for 

base  and  emitter  parasitic  resistances  that  are  critical  at 
those  currents.  Note  that  gi2  is  considered  so  small  that  in 

ordinary  SPICE  this  parameter  is  not  printed  in  the  operating 
point  output.  The  components  of  base  current  which  depend  on 
collector  emitter  voltage  are  the  recombination  in  quasi- 
neutral base  (QNB)  and  the  generation  in  collector  junction 
space-charge  region  (SCR)  [Sz81] . These  current  components 
are  typically  negligible  for  contemporary  vertical  BJTs,  and 
tend  to  reduce  base  current  with  increasing  Vce*  In  principle, 
modulation  in  QNB-width  modulates  the  base-spreading 
resistance  also,  so  that  base  current  decreases  with 
increasing  Vce*  Therefore,  gi2  predicted  by  electrical-only 

model  is  negative  and  small.  Impact  ionization  can  affect  the 
base-current  significantly  [Po89],  especially  for  advanced 
devices,  but  again,  drives  gi2  to  negative  value. 

Figure  5.5  shows  measured  and  simulated  value  of  Y22/ 
both  magnitude  and  phase  for  Ic  = 1.5  mA  and  Vce  = 6 V.  The 
discrepancies  in  the  phase  of  Y22  at  high  frequencies  (~  107 
Hz)  are  probably  the  measurement  errors  (e.g.,  components 
become  more  inductive).  Figure  5.5a  shows  that,  in  this 


Reverse  Transconductance 
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Collector  Current  Ic  (A) 


Fig.  5.4.  Common-emitter  reverse  transconductance  as  a 

function  of  collector  current,  measured  and 
simulated.  Unmodified  SPICE  gives  a small  negative 
value  for  this  parameter. 
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Fig.  5.5.  Measured  and  simulated  output  admittance.  (a) 

Magnitude,  (b)  Phase.  Vce  = 6 V and  Ic  = 1.5  mA. 
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transistor,  thermal  effects  disappear  for  frequencies  higher 
than  about  300  kHz. 

As  Figs.  5.2  - 5.5  show  the  thermally  modified  SPICE 

r*- 

simulations  agree  well  with  measurements,  demonstrating  the 
superiority  of  the  new  model  compared  to  the  electrical-only 
model.  Since  the  measurements  and  simulations  show  good 
agreements  on  the  device  level,  in  the  following  section, 
circuit  simulation  results  predicted  by  the  modified  SPICE 
will  be  assumed  to  be  the  same  as  measurements. 

5.3.  Thermal  effects  in  integrated  circuits 
The  large-signal  effects  of  self-heating  on  IC 

performance  will  be  considered  first.  Then  the  effects  on 

► 

small-signal  performances  will  be  considered.  Circuit 
simulations  will  also  be  included  for  verification  and 
demonstration.  The  NPN  transistors  used  for  the  simulations 
had  the  following  specifications:  emitter  length  of  4 (1m, 

width  of  10  |im,  collector-base  junction  depth  of  0.4  (im,  and 
Nepi  = 7xl015  cm-3.  The  value  Va  = 150  V was  extracted  at 
low  currents.  This  device  had  about  the  same  current -driving 
capabilities  as  the  23x23  (p.m)2  transistor  used  in  Section 

5.2  as  this  device  was  more  heavily  doped  than  the  other. 
However,  the  thermal  resistance  of  this  transistor  was  290 
°C/W,  higher  than  the  7 0 °C/W,  Rth  of  the  other  device. 

Therefore,  more  significant  thermal  effects  are  expected  in 
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circuits  using  this  device  than  the  one  used  in  Section  5.2. 

The  thermal  impedance  model  presented  in  Chapter  3 is 
not  directly  applicable  for  lateral  PNP  transistors.  However, 
it  is  valid  for  vertical  PNP  transistors.  As  discussed  in 
Chapter  4,  vertical  PNP  transistors  have  been  shown  to 
exhibit  similar  thermal  spreading  resistances  as  NPN  for  a 
given  emitter  area.  In  the  simulations  to  follow,  the  model 
parameters  of  the  4x10  (|i.m)2  NPN  transistor  are  also  used  for 
the  PNP  transistors.  Due  to  lower  hole  mobility,  actual  PNP 

transistors  would  show  lower  collector  current  and  larger 
transit  time  Tp  than  the  same  dimension  NPN  for  a given  base- 

emitter  voltage.  However,  using  the  same  parameters  makes 
little  difference  on  the  simulation  results. 


5.3.1.  DC  thermal  effects  in  IC 

The  increment  in  emitter- junction  temperature  over  the 
ambient  due  to  collector  power  dissipation  is  given  by  At  = 

RthIcVcE  so  that  the  DC  thermal  effects  are  increasingly 
significant  toward  higher  Rth»  Ic»  and  Vce-  Typically,  the 
terminal  currents  increase  roughly  7%/°C  when  the  base- 
emitter  voltage  is  fixed,  and  roughly  0.5  %/°C  when  the  base 
current  is  fixed  [Sz81] . For  the  current-driven  case,  the 
base-emitter  voltage  varies  roughly  -2  mV/°C.  More  exactly, 
the  temperature  coefficient  of  terminal  currents  and  base- 
emitter  voltage  is  a decreasing  linear  function  of  Vbe/  given 


/ 


by  (4.4),  (4.8),  and  (14)  in  [S176],  Therefore,  thermal 

effects  are  more  significant  at  low  Vbe  for  the  same  amount  of 

power  dissipated. 

With  regard  to  circuit  design,  the  problems  caused  by 
the  DC  thermal  effects  can  generally  be  classified  into  two 
categories.  First,  thermal  effects  magnify  the  dependence  of 
collector  current  on  Vqe/  and  second,  thermal  effects  make  the 

V' 

terminal  currents  depart  from  the  exact  exponential 

/ 

dependence  on  the  base-emitter  voltage. 

The  significance  of  the  former  effect,  which  is  most 
significant  for  the  voltage-driven  case,  will  first  be 
considered  through  the  example  of  a current  mirror.  The 

J 

current  mirror  is  an  ubiquitous  building  block  in  analog 
integrated  circuits.  Figure  5.6  shows  a basic  current-mirror 
circuit,  which  is  susceptible  to  an  offset-current  (Iout  - 
Iref)  error  not  only  due  to  the  limited  Early  voltage,  but 

also  due  to  local  thermal  feedback. 

Assuming  the  transistors  are  perfectly  matched,  an 
offset  current  can  be  developed  due  to  unequal  voltages  on 
the  collectors  of  Qi  and  Q2 . An  electrical-only  BJT  model 
would  predict  an  offset  current  error  resulting  from  the 
difference  in  base-width  modulation  while  the  actual  value  of 
offset  current  will  also  be  affected  by  thermal  effects.  For 
lref  = 1 mA  and  Vc2  - Vci  = 10  V,  the  actual  temperature  of  Q2 
is  higher  than  Qi  by  roughly  Rth^ref  (Vc2  ~ Vci)  which  is  about 
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Fig.  5.6.  Circuit  diagram  of  a simple  current  mirror. 
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2.9  °C.  This  is  effectively  the  same  condition  as  the  base- 
emitter  voltage  of  Q2  is  driven  by  roughly  5.8  mV  (Vbe  varies 
by  about  -2  mV/°C)  harder  than  Qi . Figure  5.7  shows  the 
offset  current  (lout  ~ Iref)  as  a function  of  Iref  for  a 
differential  collector  voltage  (Vc2  ~ Vci)  of  10  V based  on 
modified  and  ordinary  SPICE  simulations.  As  can  be  seen  in 
Fig.  5.7,  the  errors  between  the  two  offset  currents  become 
increasingly  significant  toward  higher  differential  power.  At 
Ic  = 1 mA,  ordinary  SPICE  predicts  an  offset  current  of  75  |1A 

while  the  actual  value  is  226  |iA,  a ratio  of  about  1 to  3 . 

In  above  example,  one  might  suggest  that  the  error  can 
be  avoided  by  extracting  Va  at  the  actual  operating  point  of 

those  transistors.  However,  the  extracted  Early  voltage 
varies  strongly  as  a function  of  Ic  due  to  thermal  effects. 

Early  voltage  can  even  change  significantly  over  the 
collector-emitter  voltage  swing.  Therefore,  extracting  Va  at 

the  actual  operating  current  is  not  practical  for  the  purpose 
of  predicting  circuit  performance. 

In  fact,  thermal  modeling  is  necessary  to  avoid  errors 
even  for  the  simpliest  case.  Fig.  5.8  show  a measured  Ic 
versus  Vce  of  the  4x10  (|im)  NPN  transistor,  for  the  cases  base 
is  driven  by  a voltage  source  and  current  source.  As  can  be 
seen  in  the  figure,  collector  currents  which  are  same  at  low 

Vce  diverge  significantly  toward  high  Vce-  These  diversions 

/ # 

are  due  to  the  discrepancies  in  temperature  coefficient  of 


Offset  Current  Iout  - Iref  (mA) 
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Reference  Current  Iref  (mA) 


Fig.  5.7.  Offset  current  (lout  ~ Iref)  versus  mirror  current 

(Iref)  of  the  simple  current  mirror  for 
differential  collector  voltages  (Vc2  - Vci)  of  10 
V. 


Collector  Current  I,-,  (mA) 
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Fig.  5. 


Collector-Emitter  Voltage  VCE  (V) 


8 


Typical  measured  collector  currents  as  a function 
of  collector-emitter  voltage.  Emitter  area  = 4x10 

(Jim)2  and  Nepi  = 7xl015  cm-3. 
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collector  current  which  depend  on  the  bias  conditions. 
Therefore,  even  at  a given  collector  current  Ic,  a single 
value  Va  can  not  satisfy  the  output  characteristic  as  the 
small-signal  slope  of  the  output  depends  on  the  bias 
condtion. 

Fig.  5.9  shows  a self-biasing  VBE_based  current 
reference  which  is  another  circuit  susceptible  to  the  Early 
voltage  degradation.  In  Fig.  5.9,  changing  the  supply 
voltage,  Vcc>  mostly  varies  the  collector-base  voltage  of 
transistor  Q3 , Vce3 # and  thermal  effects  make  the  current 

mirror  output  Ic3  to  be  strongly  dependent  on  Vce3  • For  vcc  = 

/ 

12  V,  ordinary  SPICE  predicts  0.09  (0.003)  for  the 
sensitivity  of  the  output  current  lout  to  the  supply  voltage 
Vcc  which  is  defined  as  (Vcc/lout)  (dlout/^Vcc)  • For  t-*ie  same 
circuit,  modified  SPICE  predicts  0.45  (0.021),  which  is  more 
than  five  times  higher  than  what  the  ordinary  SPICE  predicts. 
This  amount  of  error  can  turn  a useful  circuit  into  a totally 

] 

inadequate  one. 

In  both  of  above  examples,  increasing  the  base  or 
emitter  resistance  can  help  reduce  the  errors.  The  negative 
electrical  feedback  across  Rb  or  Re  opposes  the  increase  in 
terminal  currents.  This  will  be  explained  more  detail  in  the 
section  on  small-signal  analysis.  In  the  supply- independent 
current  source,  replacing  the  current  mirror  pair  Q3  and  Q4 
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Fig.  5.9.  Circuit  diagram  of  a self-biasing  Vbe  reference. 
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with  a cascode  or  Wilson  current  source  can  greatly  reduce 

the  sensitivities.  In  a cascode  source,  such  as  in  Fig.  5.10, 

* 

the  output  current  is  mostly  determined  by  the  common-emitter 
transistor  Q3 . Changing  the  output  voltage  Vc2  primarily 
changes  the  power  in  common-base  transistor  Q2 . Since  Q2  is 

driven  with  nearly  constant  current,  the  self -heating  of  Q2 

# *• 

mostly  changes  its  base-emitter  voltage,  and  has  little 
effect  on  the  output  current. 

A class  of  circuits  susceptible  to  the  deviations  from 
the  exponential  dependence  of  collector  current  to  the  base- 
emitter  voltage  are  the  translinear  circuits.  Translinear 
circuits  exploit  the  near-exact  proportionality  of  the 
transconductance  of  a BJT  to  its  collector  current  (or  the 
near -exponential  dependence  of  the  collector  current  on  base- 
emitter  voltage)  hence  the  term  translinear. 

Figure  5.11  shows  a four-quadrant  multiplier  cell. 
Depending  on  the  input  voltages  vi  and  V2,  collector  currents 
will  vary  between  0 and  2 mA.  Therefore,  the  temperature  of 

m 

each  transistor  can  vary  substantially  as  a function  of 
inputs.  The  variations  in  the  temperature  of  each  device 
introduce  extra  non-proportionality  of  transconductance  to 
its  collector  currents  so  that  the  differential  output 
current  (Ii  - I2)  deviates  more  from  linear  behavior.  Figure 

5.12  compares  the  deviation  of  differential  output  current 
(II  - I2)  from  the  linear  behavior  as  a function  of  input 
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Fig.  5.10  Circuit  diagram  of  a cascode  current  source. 
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Fig.  5.11  Circuit  diagram  of  a four  quadrant  multiplier. 


Linear  Function  (mA) 
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Input  v-l  (V) 


Fig.  5.12.  Deviation  of  the  differential  output  current  (Ii  - 

I2)  from  the  linear  behavior  as  a function  of  input 
vi  for  V2  = 0.5  V. 
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voltage  vi  for  V2  = 0.5  V,  based  on  ordinary  and  modified 
SPICE  simulations.  As  can  be  seen  in  Fig.  5.12,  thermal 
effects  make  the  differential  output  current  a nonlinear 
function  of  the  inputs  in  the  region  where  ordinary  SPICE 
predicts  a linear  behavior.  This  is  direct  indication  of  a 
large  amount  of  distortion. 

Figure  5.13  shows  a typical  band-gap  voltage  reference 
circuit.  In  the  figure,  the  emitter  area  of  Qi  is  ten  times 
that  of  Q2 . Resistor  R2  is  made  to  be  10-R3,  so  that  Ic2  = 
10-Ici-  Note  that  the  ten  times  smaller  transistor  Q2 
dissipates  ten  times  more  power  than  Qi . If  the  two 
transistors  Qi  and  Q2  are  assumed  to  operate  at  the  same 
temperature,  the  output  voltage  Vref  can  be  expressed  as 


Vref  = Vbei  + Vt  [ln( 


Ie2  Isi 

iEl  IS2 


) (1  + 


(5.2) 


The  value  of  Vref  is  significant,  as  it  is  the  point  where  the 
negative  and  positive  temperature  coefficient  of  the  first 
and  second  term  in  (5.2)  cancels  out,  so  that  the  temperature 
coefficient  of  reference  voltage  becomes  zero  at  Vref.  This 

circuit  is  designed  to  have  zero  temperature  coefficient  at 
temperature  T = 0 °C  with  Vref  = 1.293  V.  Figure  5.14  shows 

the  simulation  results  based  on  ordinary  (dashed  line)  and 
modified  (solid  line)  SPICE  for  Ic2  — 1 mA  and  Vcc  = 5 V. 

Thermal  effects  lead  to  about  10  mV  discrepancies  in  the 
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Fig.  5.13.  Circuit  diagram  of  a band-gap  voltage  reference. 
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Temperature  (°C) 


Fig.  5.14.  Output  voltage  of  band-gap  voltage  reference  as  a 

function  of  temperature . 
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output  voltage  Vref  and  5 °C  displacement  in  T0,  which  is 
defined  as  the  temperature  at  which  dvref/3T  = 0.  These  errors 
are  caused  by  the  deviation  of  Ic  from  the  near-exponential 
VBe  dependence.  A popular  way  to  tune  the  temperature 

coefficient  of  the  band-gap  reference  is  to  trim  the  output 
to  the  value  of  Vref  [Ku73,  Br74],  assuming  the  right  output 

value  will  give  the  right  temperature  coefficient.  Here, 
resistor  Ri  is  adjusted  to  match  the  output  voltage  at  T = 0 
°C  to  the  value  predicted  by  ordinary  SPICE  (1.293V).  Figure 
5.14  (dotted  line)  shows  the  result  of  modified  SPICE 
simulation  after  Ri  has  been  adjusted.  Corrections  without 
understanding  the  error  source  can  worsen  the  errors.  As  can 
be  seen  in  Fig.  5.12,  this  procedure  actually  displaces  T0  by 

more  than  30  °C . 

In  this  section,  circuit  implication  of  DC  thermal 
effects  have  been  studied.  The  behavior  of  device  level 
thermal  effect  that  can  give  significant  errors  in  many 
analog  circuits  are  the  strong  dependence  of  Ic  on  Vce  and  Ic 
being  a less-exact  exponential  function  of  Vbe,»  These 
behaviors  lead  to  errors  in  many  precision  analog  circuits 
such  as  increasing  the  offset  currents,  degrading  the  supply 
voltage  sensitivity,  distorting  the  output  behavior,  and 
degrading  the  temperature  dependence  of  output  voltage.  In 
general,  large-signal  thermal  effects  are  insignificant  if 
the  power  dissipation  is  small.  The  errors  caused  by  stronger 
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dependence  of  Ic  on  Vce  can  be  reduced  considerably  through 
electrical  negative  feedback  or  cascoding. 

I 


5.3.2.  Small-signal  thermal  effects  in  IC 

4 

In  small-signal  analysis,  the  thermally  induced 
variations  in  terminal  currents  or  voltages  lead  to  change  in 
the  y-parameters . The  common-emitter  y-parameters  corrected 
for  thermal  feedback,  (3.6),  are  restated  here 


mn 


YmnE  + Dm^thlmln 
1 - DmZthP 


(5.3) 


where  m or  n = 1 for  the  base  or  2 for  the  collector,  YmnE  is 
the  uncorrected  electrical  Y parameter,  Zth  is  the  thermal 
impedance,  and  P is  the  dissipated  power.  D2  is  the 

fractional  temperature  coefficient  of  the  collector  current, 
and  Di  is  that  of  the  base  current.  Typically,  at  low 
frequencies  Y12E  and  Y22E  are  both  small  so  that  the  thermal 
term,  DmZthImIn»  can  easily  dominate.  Conversely,  Yue  and  Y21E 
are  larger  to  start  with,  so  that  thermal  effects  on  these 
parameters  are  significant  only  for  high  powers  in  which  case 
the  denominators  become  considerably  less  than  one.  For 
example,  for  the  4x10  (|im2)  transistor  used  here,  the 
denominator  of  (5.3)  becomes  10%  less  than  1 at  Ic  = 1.5  mA 

assuming  D2  = 4.5  %/°C  and  Vce  = 5 V.  At  Ic  = 


1.5  mA,  the 
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second  term  of  (5.3)  is  only  about  6xl0-4%  of  the  first  term 
for  both  gn  and  g2i- 

The  denominator  of  (5.3)  is  actually  a measure  of  the 
significance  of  DC  thermal  effects.  DmRthP  term  includes  all 
the  parameters  determining  the  DC  thermal  effects.  If  DroRthP 
approaches  1,  the  transistor  will  get  into  thermal  runaway. 
Thermal  runaway  requires  very  high  powers  and  typical 
integrated  circuits  are  unlikely  to  operate  in  that  region, 
unless  the  thermal  impedances  are  extremely  high.  Thermal 
runaway  is  an  important  aspect  of  thermal  effects,  and  has 
been  well  known  in  power  electronics  [Sc67,  Ho74,  Ar74, 
La81] . This  work  excludes  such  a high  power  local-heating 
effects.  Therefore,  for  the  modest  power,  assuming  (l-ImRthp) 

= 1,  the  significance  of  small-signal  thermal  effects  is 

jj  ^ 

determined  by  the  size  of  Rth»  Ib/  Ic/  and  the  temperature 
coefficient  of  Ib  and  Ig- 

Thermal  effects  on  the  low-frequency  small-signal 
behavior  of  integrated  circuits,  will  be  considered  here 
first.  Self-heating  affects  y-parameters  the  most  at 
frequencies  close  to  DC,  where  (5.3)  simplifies  to 

9nm  — 9mnE  + PmRth^mln.  (5.4) 

As  mentioned  above,  gi2  and  g22  are  the  DC  small-signal 
parameters  most  affected  by  thermal  feedback.  As  mentioned 
previously,  any  other  physical  mechanisms  except  thermal 
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effects  tend  to  make  gi2  negative.  However,  thermal  effects 
make  gi2  positive  even  at  low  to  medium  currents  since  gi2E  is 

so  small  in  today's  typical  BJTs. 

Figures  5.15a  and  5.15b  show  a one-transistor  common- 
emitter  voltage-gain  stage  and  its  hybrid-Jt  small-signal 

equivalent  circuit.  In  Fig.  5.15b,  r^  = 1/gn,  gjj.  = gi2/  Om  = 
g2 1 , and  rQ  = l/g22,  so  that  thermal  effects  can  be 

incorporated  into  these  conductances  based  on  (5.4),  and  Re, 
Rb,  and  Rc  represent  the  total  resistance  looking  out  at 

emitter,  base,  and  collector  port,  respectively,  including 
parasitic  resistances.  From  (5.4),  rG  can  be  expressed  as 

roE (c  c] 

y ° ~ 1 + D2RthIcVA 

where  r0E  is  the  electrical  output  resistance  which  is  defined 
as  VA/ Ic  • Here  VA  represents  the  Early  voltage  measured  at  low 
currents  where  the  output  conductance  is  controlled  entirely 
by  electrical  effects.  Note  that  the  dimensionless  term 
D2RthIcVA  in  (5.5)  is  a measure  of  the  significance  of  thermal 
effects  on  rQ.  This  term  involves  all  the  factors  determining 
the  thermal  degradation  of  rD.  The  greater  this  term,  the 
more  rQ  is  degraded.  rQ  drops  to  half  of  what  the  electrical- 
only  model  predicts  when  D2RthlcVA  = 1. 

From  Fig.  5.15b,  the  output  resistance  at  the  collector 


is  given  by 
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Fig.  5.15.  A common-emitter  voltage  gain  stage,  (a)  Circuit 

diagram,  (b)  Equivalent  circuit  diagram. 
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Rout  = rQ  + Rg/  /Rb  + 


] (; 


^7tRE 


Re  + gfxro(RE  + Rb) 

tcji  + (Re  + Rb)  (1  - Po9|a.ro)  + 


Po^o) 


(5.6) 


where  po  is  assumed  to  be  equal  to  gmrji.  Based  on  (5.4), 

assuming  gi2E  is  negligible,  and  assuming  the  small-signal 
current  gain  pQ  equals  the  DC  gain  Iq/Ib / the  l-por0g|j.  term  in 

(5.6)  can  be  expressed  as 


i o _ 1 + DBVARthic 

1 P0r0gH  - i + D2VARthIc 


(5.7) 


where  Dp  is  the  temperature  coefficient  of  (3.  Typically, 

DpVARthlc  « D2VARthIc  as  Dp  « D2 . As  can  be  seen  from  (5.7), 
the  (1  - p0r0g|x)  term  is  always  greater  than  zero. 

From  (5.6),  for  RE  and  RB  « rn,  the  difference  between 
Rout  and  the  output  resistance  predicted  by  electrical-only 
model,  R0utE-  is  controlled  by  the  degradation  in  rQ.  In  this 
case,  measured  output  resistance  can  be  significantly  lower 
than  ordinary  SPICE  predicts  depending  on  the  size  of 
D2RthIcVA.  With  increasing  RB  and  RE,  the  error  tends  to  be 
reduced.  For  RB  or  RE  » rn,  Rout  stays  about  the  same  as 
RoutE;  Rout  drops  only  due  to  the  variation  in  P with 

temperature  which  can  be  considered  as  a second-order  thermal 
effect.  For  example,  for  the  limiting  case  of  Rb  — » 0 and  Re 

» rK  (i.e.,  near-constant  emitter  current),  (5.6)  simplifies 
to 

1 + Po  _ RoutE 
1 - Por0g^  ~ 1 + DpRthIcVA 


Rout  — rQ 


(5.8) 
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where  R0utE  = roE(l+P0)*  As  discussed  above,  typically,  the 
DRRthlcvA  term  is  substantially  smaller  than  1.  In  (5.8),  the 
degradation  in  rQ  is  canceled  out  to  first  order  by  the  (1  - 
Po^ogfi)  term. 

The  input  resistance  Rin  of  the  same  circuit  is  given 


Rin  — r7t  + Re/ /Rc  + 


Re  + Po  (Rf,  + Rc) 


_rG  + (Re  + Rc ) ( 1 - Po  r0) 


r0RE 

,Re  + Rc 


+ rjtrog^ 


(5.9) 

Note  that  this  expression  can  be  obtained  from  (5.6)  by 
inspection,  by  exploiting  the  symmetry  of  the  equivalent 
circuit.  Notable  in  (5.9)  is  the  strong  dependence  of  Rin  on 

the  loading  at  the  output  port.  For  example,  for  the  limiting 
case  of  Re  ^ 0,  (5.9)  simplifies  to 


rn 

P0g|x(r0//Rc) 


(5.10) 


For  Rc  « rD,  (5.10)  simplifies  to  Rin  = rrt.  Thermal  effects 
on  rjj  are  negligible.  However,  in  the  other  limit  Rc  >>  rQ, 
the  denominator  of  (5.10)  can  be  substantially  smaller  than 
1,  so  that  the  input  resistance  can  be  many  times  greater 
than  what  the  electrical-only  model  predicts.  This  effect  can 
be  explained  as  follows:  an  increment  in  the  input  voltage 
tends  to  increase  the  collector  current;  this  incremental 
current  Ale  reduces  the  collector-emitter  bias.  For  large 
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enough  Rc  (i.e.,  Rc  » rQ)  , the  power  dissipation  can  be 
reduced  even  though  the  collector  current  has  been  increased. 
In  that  case,  the  base  current  increases  less  than  that 
predicted  by  the  electrical-only  model,  as  the  emitter- 
junction  temperature  actually  drops,  increasing  the  effective 
input  resistance. 

These  thermally  induced  variations  in  input  and  output 
resistances  lead  to  errors  in  other  circuit  responses  such  as 
voltage  ’ The  lower  output  resistance  can  significantly 


reduce  the  low-frequency  voltage  gain.  In  some  cases,  the 
increase  in  input  impedance  can  partially  cancel  the  effect 
of  decreasing  output  impedance.  In  such  a case,  the  net 
effect  on  gain  can  be  small  even  when  both  impedances  are 
strongly  in  error.  Figure  5.16  shows  an  equivalent  circuit  of 
such  case.  In  Fig.  5.16,  Rs  is  the  output  resistance  of 
previous  stage,  Rin  the  input  resistance,  gm  the 
transconductance.  Rout  the  output  resistance,  and  Rl  the  load 
resistance  of  amplifying  stage.  Under  the  conditions  that  Rs 
» Rin  and  Rl  » Rout  (that  is,  the  gain  stage  is  driven  by  a 
high  impedance  source  and  has  a high  impedance  load) , the 
voltage  gain  of  this  circuit  is  given  by 


(5.11) 
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Fig.  5.16.  Equivalent  circuit 

stage. 


diagram  of  a typical  gain 
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In  (5.11),  the  higher  Rin  and  the  lower  Rout  are  multiplied 
together  so  that  the  net  value  stays  about  the  same  as  what 
the  electrical-only  model  predicts.  This  is  an  important  case 
as  such  situation  happens  in  the  gain  stage  of  typical  two- 
stage  op-amps  (e.g.,  741).  Even  though  there  is  little  effect 
on  the  DC  gain,  errors  in  the  input  and  output  impedance  can 
lead  to  significant  errors  in  the  frequency  response  of  the 
circuit,  as  will  be  discussed  later. 

The  small-signal  aspects  of  thermal  effects  so  far  are 
demonstrated  in  following  example.  Figure  5.17  shows  a 
common-emitter  amplifier  with  a PNP  active  load.  By 
inspection,  the  output  resistance  of  this  circuit  is 
determined  by  the  parallel  combination  of  the  output 
resistances  of  transistor  Qi  and  Q2 . From  the  discussions 

above,  the  gain  of  this  circuit  can  be  considerably  less  than 

what  electrical-only  model  predicts.  Figure  5.18  shows  the  DC 

transfer  characteristic  of  this  circuit  based  on  modified  and 
ordinary  SPICE  simulations  for  lCi  = 1 mA.  The  figure  shows 

that  thermal  feedback  not  only  reduces  the  DC  voltage  gain 
but  makes  it  a more  strongly  nonlinear  function  of  the  input 
swing.  This  is  caused  by  the  strong  Vce  dependence  of  the 
current -mirror  output  Ic2  • From  the  modified  SPICE  simulation, 
Ic2  varies  about  50%  over  the  output  swing. 

To  confirm  the  thermal  effect  on  the  input  resistance, 
consider  the  circuit  of  Fig.  5.17  with  the  simple  current 
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Fig.  5.17.  Circuit  diagram  of  a simple  common- emitter  with 

active-load. 


Output  Vout 
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Input  Vin  (V) 


Fig.  5.18.  Transfer  characteristic  of  a simple  common -emit ter 

circuit . 
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mirror  (Q2  - Q3)  replaced  with  a Wilson  current  source  (Fig. 
5.19).  In  this  case  the  output  resistance  of  the  circuit  is 
given  by  (5.10).  Modified  SPICE  predicted  an  input  resistance 
of  1.5xl04  Q while  ordinary  SPICE  gives  6.1xl03  Q,  a 
prediction  error  of  146  %. 

Figures  5.20a  and  b show  the  gain  and  phase  responses  of 
the  circuit  in  Fig.  5.17.  In  both  gain  and  phase,  the  results 
of  simulations  with  and  without  thermal  effects  show 
considerable  discrepancies.  The  low-frequency  voltage  gain 
predicted  by  modified  SPICE  is  1.49xl03  while  ordinary  SPICE 
predicts  4.53xl03,  a prediction  error  of  over  200  %. 

Note  that,  in  Fig.  5.20a,  the  voltage  gain  from  modified 
SPICE  varies  at  frequencies  where  ordinary  SPICE  predicts  a 
constant  value,  especially  near  lxlO6  Hz.  As  can  be  seen  in 
(5.3),  the  frequency  dependence  of  Zth  in  y-parameters  can 

cause  the  input  impedance,  output  impedance,  and  the  gain  of 
a circuit  to  vary  with  frequency  at  frequencies  where 
ordinary  SPICE  predicts  a constant  value.  The  frequency 
response  of  Zth  depends  on  the  device  structure  [Fo91b] . The 
narrower  the  base-width  and  the  emitter  area,  the  higher  the 
frequency  over  which  thermal  feedback  will  present.  From 
Chapter  3,  in  a uniform  medium,  the  time  at  which  the  impulse 
response  at  a point  a distance  r away  reaches  maximum  is 
given  as  tpk  = r2/6K.  At  tpk  the  temperature  increases  to 

roughly  10%  of  its  final  value.  Therefore,  in  principle,  the 
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Fig.  5.19  Circuit  diagram  of  a simple  common-emitter  with  a 

Wilson  source  as  active-load. 


Phase  ( ) Voltage  Gain 
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Frequency  (Hz) 
(a) 


Frequency  (Hz) 
(b) 


Fig.  5.20.  Frequency  response  of  a simple  common-emitter 

circuit  with  active-load,  (a)  Gain,  (b)  Phase. 
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infinitesimal  heat  source  located  in  collector- junction  SCR 
will  raise  the  temperature  of  the  closest  emitter  junction  to 
10%  of  its  final  value  within  Wb2/6k  seconds  (Wb  is 

metallurgical  base-width) , while  for  the  farthest  emitter 
junction,  it  takes  [ (Wb+H) 2+W2+L2 ] / 6K  seconds.  Typically, 

base-width  is  much  smaller  than  emitter-width  and  length  so 
that  the  emitter  area  is  the  dominating  factor  to  determine 
tpk  of  device,  or  the  frequency  characteristic  of  Zth-  From 
(5.1b)  and  (5.1c),  reff  = 2VwLfif2- 

In  Fig.  5.20b,  the  hump  in  the  phase  response  of  the 
modified  SPICE  simulation  is  due  to  the  negative  phase  of  Zth 
in  Y22  as  seen  in  Fig.  5.5b. 

The  stability  of  a circuit  can  be  affected  by  the 

\ : I • ‘"y  y •'  W 

frequency  dependence  of  the  thermal  impedance  discussed 

above.  At  low  frequencies,  the  stability  of  a circuit  can 
also  be  affected  more  significantly  by  the  thermally  induced 
modulations  in  input  and  output  resistances.  The  effect  of 
the  modulation  in  input  and  output  resistance  to  the  pole 
frequencies  of  a circuit  will  be  considered  here.  Similar 
arguments  can  be  applied  to  zeros  also.  Assume  that  the  pole 
frequencies  are  low  enough  that  (5.4)  is  valid  and  the 
collector  current  is  high  enough  to  cause  significant  thermal 
effects.  Depending  on  the  circuit  configuration,  four  cases 
are  possible.  If  pole  frequency  fp  is  controlled  by  an  output 
resistance,  the  smaller  R0ut  due  to  the  thermal  effects  can 
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increase  fp  compared  to  the  electrical-only  model  prediction. 
If  this  fp  was  the  dominant  pole  frequency,  the  circuit  will 
have  smaller  phase  margin  than  ordinary  SPICE  predicts, 
assuming  the  rest  of  the  poles  are  not  affected  by  thermal 
feedback  or  by  the  variation  in  this  pole.  The  opposite  can 
be  true  if  this  fp  was  the  frequency  of  a nondominant  pole. 
For  the  case  where  the  pole  frequency  fp  is  controlled  by  an 
input  resistance,  which  can  be  raised  by  self -heating,  the 
thermal  feedback  moves  the  pole  to  lower  frequency.  Again, 
depending  on  whether  the  pole  was  dominant  or  not,  this 
variation  in  fp  can  increase  or  decrease  the  phase  margin  of 

the  circuit.  Table  1 summarizes  the  resulting  effects  for  the 
four  possible  cases. 

A change  in  pole  frequency  by  less  than  an  order  of 
magnitude  has  little  effect  on  the  phase  margin  of  the 
circuit  if  the  rest  of  the  poles  were  initially  separated  by 
more  than  two  decades.  At  frequencies  more  than  a decade  from 
the  pole  frequency,  phase  effects  are  small.  For  Cases  (a) 
and  (d)  of  Table  1,  if  the  pole  influenced  by  thermal 
feedback  was  initially  separated  by  less  than  two  decades 
from  some  other  pole,  thermal  feedback  can  degrade  the  phase 
margin.  In  Cases  (b)  and  (c)  , the  circuit  will  be  over- 
compensated. Considering  the  chip  area  that  will  be  taken  up 
by  the  compensation  capacitors,  accurate  prediction  of 
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Table  1:  List  of  various  conditions  under  which  thermal 

feedback  can  affect  the  stability  of  a circuit 
resulting  from  the  discrepancies  in  the  pole 
frequency  fp. 


Effect  on  fD 

If  pole  constitutes 

Resulting  effects 

Smaller  R0ut 
increases  fp 

(a)  Dominant  pole 

Reduces  phase  margin 

(b)  Nondominant  pole 

Increases  phase  margin 

Larger  Rin 
reduces  fp 

(c)  Dominant  pole 

Increases  phase  margin 

(d)  Nondominant  pole 

Reduces  phase  margin 
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frequency  response  is  required  to  optimally  size  the 
compensating  capacitor. 

Note,  as  discussed  in  the  DC  aspects,  correcting  the 
forward  Early  voltage  Va  cannot  solve  the  problems  induced  by 
thermal  effects.  If  Va  is  extracted  at  low  current, 
simulation  results  are  accurate  at  low  current  for  all 
frequencies  and  at  high  frequencies  for  all  currents. 
However,  simulation  results  will  be  wrong  at  low  frequencies 
with  moderate  to  high  currents.  On  the  other  hand,  if  Va  is 

extracted  at  high  currents,  simulation  results  will  be 
accurate  at  high  currents  for  low  frequencies  but  wrong  at 
low  currents  for  all  frequencies  and  at  high  frequencies  for 
all  currents.  Either  way,  designers  cannot  avoid  errors 
unless  thermal  effects  are  included  properly  in  the 
simulator.  Table  2 summarizes  the  above  discussions. 

A possible  mistake  in  circuit  compensation  by  not 
understanding  the  aspect  summarized  in  Table  2b  is 
illustrated  in  the  following  example.  Figure  5.21  shows  a 
two-stage  amplifier  in  which  the  second  stage  will  generate  a 
pole  at 

IPil  = (5*12) 

•T'-pCp 

where  Rp  and  Cp  are  the  capacitance  and  resistance  of  the 
second  stage,  respectively.  For  Rc  » rQ,  based  on  the  BJT 
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Table  2.  Summary  of  the  relations  between  ordinary  SPICE 
simulations  error  and  operating  conditions  depending  on 
whether  the  Early  voltage  is  extracted  at  high  currents  or 
low  currents.  Extracting  Va  at  low  currents  represents  Va 
considering  only  the  electrical  effects,  high  frequencies 
represents  high  enough  frequencies  to  avoid  thermal  effects, 
and  low  currents  represents  negligible  power  dissipation. 


(a)  Va  extracted  at  low  Ic 


Simulation  results 

Low  currents 

High  currents 

Low  frequencies 

Accurate 

Wrong 

High  frequencies 

Accurate 

Accurate 

(b) Va  extracted  at  high  Ip 


Simulation  results 

Low  currents 

High  currents 

Low  frequencies 

Wrong 

Accurate 

High  frequencies 

Wrong 

Wrong 
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Fig.  5.21.  Small-signal  equivalent  circuit  diagram  of 

stage  amplifier. 


two 
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equivalent  circuit  shown  in  Fig.  5.15b,  the  first  stage  has  a 
pole  near 


Cji  (Rb  / / rrc)  + c cs 


i 

rQ  + gmr0(RB//rJ 

11  1 - gmr0gn(RB//r7i) 


(5.13) 


where  cn,  c^,  and  ccs  are  the  parasitic  capacitances  of 
transistor  Qi.  Let's  assume  I pi  I < IP2U  and  I p2  I is  high 
enough  that  thermal  effects  are  negligible  at  that  frequency, 
and  assume  this  circuit  need  to  be  compensated  for  unity  gain 
with  a phase-margin  of  45°.  One  could  add  a capacitance  Cp'  in 
parallel  with  Cp  to  lower  the  dominant  pole  frequency.  Cp' 
must  satisfy  the  relation 


Ayp 

Rp  I P2  I 


GmlRplGp 


IP2  I 


(5.14) 


where  Avo  is  the  low-frequency  voltage  gain  of  the  circuit, 
Gmi  the  transconductance  of  first  stage,  R0i  the  output 
resistance  of  first  stage.  I P2  I and  R0i  must  be  known  to 
estimate  the  size  of  Cp‘  and  designers  typically  use 
simulation.  From  (5.13),  I p2 1 is  a function  of  rD,  which  is 
directly  proportional  to  Va,  and  R0i  = r0.  Suppose  Va  of  Qi 
was  extracted  at  the  actual  operating  current  so  that  the 
low-frequency  gain  is  predicted  correctly.  Figure  5.22a  and  b 
compare  the  frequency  response  of  this  circuit  based  on  both 
ordinary  and  modified  SPICE  simulation.  In  Fig.  5.22,  Va  = 56 
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Fig.  5.22.  Gain  and  phase  versus  frequency  for  a 

amplifier. 


two  stage 


Phase 
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V,  which  was  extracted  at  Ic  = 1 mA,  was  used  for  the 

ordinary  SPICE  simulation  and  Cp ' was  adjusted  to  give  45° 

phase-margin  at  unity-gain  based  on  (5.13)  and  (5.14).  The 
same  Cp'  was  used  for  modified  SPICE  simulation,  but  Va  was 

left  to  be  150  V,  which  is  the  Early  voltage  considering 
electrical  effects  only,  as  in  real  circuit.  (as  can  be  seen 
in  Fig.  5.22a,  ordinary  SPICE  gives  the  same  voltage  gain  as 
what  modified  SPICE  does  at  low  frequencies.  However,  with 
increasing  frequency,  the  actual  voltage  gain  rises  above  the 
prediction  of  ordinary  SPICE  as  thermal  effects  diminish.  At 
frequencies  high  enough  to  neglect  thermal  effects,  output 
resistance  rQ  is  controlled  by  the  electrical  effects  only. 
Therefore,  using  Va  = 56  V leads  to  underestimation  of  Cp‘. 
This  circuit  designed  to  have  45°  ended  up  getting  only  about 
25°  phase  margin. 


5.4.  Summary 

Depending  on  the  bias  conditions,  thermal  feedback 
increases  the  terminal  currents  (Ie,  Ib»  and  Ic)  or  reduces 
the  base-emitter  voltage  of  BJTs.  These  DC  thermal  effects 
can  lead  to  significant  errors  in  analog  circuits  as  they 
make  the  terminal  currents  more  strongly  dependant  on 
collector-emitter  voltage  and  to  be  a less-exact  exponential 
function  of  base-emitter  voltage.  Many  high-accuracy  analog 
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circuits,  especially  translinear  circuits,  are  susceptible  to 
these  DC  thermal  effects. 

In  low-frequency  analysis,  thermal  feedback 
significantly  increases  the  conductances  gi2  and  g22-  Increase 
in  g22  leads  to  significant  reduction  in  output  resistance  of 
BJT  circuit  for  the  case  transistor  is  driven  by  a low 
impedance  source,  while  the  increase  in  gi2  increases  the 
input  resistance  of  the  BJT  circuit  for  Rc  >>  rQ.  Thermal 
effects  on  input  and  output  resistances  can  lead  to  errors  in 
gain.  In  addition,  the  frequency  dependence  of  thermal 
impedance  gives  an  additional  frequency  dependence  in  input 
impedance,  output  impedances,  and  gain,  which  is  not 
predicted  by  the  existing  small-signal  model.  Discrepancies 
in  input  and  output  impedances,  both  at  low  and  high 
frequency,  can  lead  to  errors  in  stability  prediction  of  the 
integrated  circuit. 

Reducing  the  maximum  collector-emitter  voltage  of  the 
transistors  by  employing  smaller  supply  voltages  or  cascoding 
can  reduce  thermal  effects  on  the  parameters  that  are 
affected  by  the  power  dissipation.  However,  Y12  and  Y22  can  be 
significantly  affected  independent  of  power  dissipation. 
^Thermal  effects  on  terminal  currents  and  output  resistance 
tend  to  be  less  significant  with  increasing  base  or  emitter 
resistance  as  the  negative  electrical  feedback  across  these 
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resistors  partially  offsets  the  positive  thermal  feedback 
effects . ; 

Note  that  as  discussed  in  previous  chapters,  errors  such 
as  those  in  the  foregoing  examples  will  be  observed  at  much 
lower  currents  in  dielectrically  isolated  (DI) , SOI,  or 
trench-isolated  transistors.  For  example,  the  DI  BJTs  had 
about  four  times  higher  thermal  resistances  than  the  silicon 
junction  isolated  transistors  [Fo91b] . Therefore,  for  DI 
devices,  the  errors  presented  in  Section  5.3  can  be  seen  at 
Ic  less  than  250  |1A.  The  predictive  thermal  impedance  model 

presented  in  Chapter  3 is  not  directly  valid  for  DI,  SOI,  and 
trench  isolated  devices.  This  model  was  developed  for 
homogeneous  materials.  However,  since  the  DC  simulations 
solely  depend  on  the  magnitude  of  thermal  resistance, 
modified  SPICE  can  be  used  for  the  DC  simulations  of  any 
structure  as  long  as  the  thermal  resistance  is  available. 


CHAPTER  6 

SUMMARY  AND  SUGGESTIONS  FOR  FUTURE  WORK 


In  this  dissertation,  predictive  modeling  of  high- 
current  output  resistance,  and  thermal  effects  on  modeling. 
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characterization,  and  circuit  design,  which  are  increasingly 
significant  with  device  scaling  and  isolation  technologies 
using  dielectric  material,  have  been  discussed.  A study  of 
the  effects  of  carrier-velocity  saturation  on  the  output 
resistance  characteristic  of  BJTs  was  done.  Physical  analysis 
demonstrates  that  BJT  output  resistance  can  be  raised 
significantly  by  carrier  velocity-saturation  especially  at 
high  currents.  Analytical  expressions  for  normalized  BJT 
output  resistance  helps  to  understand  the  effects 
qualitatively.  On  the  thermal  part,  this  work  demonstrated 
the  increasing  significance  of  local-heating  effects  with 
device  scaling  and  identified  the  roles  of  thermal  spreading 
impedance  and  chip/package  impedance  on  BJT  thermal  behavior 
depending  on  the  power  dissipation.  The  physical  analysis  of 
BJT  behavior  corrected  for  temperature  variation  demonstrates 
the  significance  of  thermal  effects  on  BJT  small-signal 
parameters,  which  justifies  the  use  of  frequency-domain 


thermal  model. \ The  thermal  impedance  model  is  simple,  yet 
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predictive,  and  can  be  used  for  compact  modeling  in  SPICE  to 
improve  the  accuracy  of  BJT  circuit  simulation.  The 
fundamental  limits  with  standard  thermal  parameter  extraction 
methods  are  pointed  out,  and  new  methods  are  introduced  which 
are  compatible  for  any  type  of  device,  including  advanced  BJT 
structures.  From  the  study  of  circuit  implications  of  self- 
heating,  this  work  provided  general  insight  into  thermal 
effects  on  integrated  circuit  design,  identified  the  aspects 
of  thermal  effects  that  can  lead  to  significant  errors  in  a 
certain  circuit  configuration,  and  provided  design  principles 
to  reduce  thermal  effects  in  integrated  circuit  design. 

For  future  work,  the  first  suggestion  goes  to  the 
exploration  of  SPICE  BJT  model  for  temperature  dependences  of 
the  SPICE  model  parameters.  Increase  in  device  temperature 
tends  to  be  really  significant  at  currents  where  the 
transistor  goes  into  high  level  injection,  and  the  device 
performance  strongly  depends  on  the  temperature  dependences 
of  the  terminal  currents  and  voltages.  This  work  relied  on 
the  temperature  models  used  in  SPICE  to  predict  the  thermal 
effects  on  high  currents.  However,  the  accuracy  of  these 
models  has  not  been  fully  explored.  Considering  the 
empiricism  prevailing  in  the  SPICE  BJT  model,  it  is  most 
likely  that  more  physical  models  for  temperature  coefficients 
are  required  for  better  prediction  of  high-current  thermal 
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Transient  thermal  effects  have  not  been  discussed  much 
in  this  work.  One  of  the  reasons  was  that,  based  on  Denison 
and  Walter's  work  [De89],  local -heating  was  expected  to  have 
little  effect  on  digital  circuit  performance.  Another  reason 
was  that  a transient  thermal  circuit  simulator  was  not 
available.  Thermal  effects  can  significantly  degrade  the 

v ■'■■■■..  • — — 

settling  time  of  a typical  output  buffer  with  large  voltage 
pulse  applied  at  the  input,  or  a switching  stage  where  the 
voltages  are  switched  by  activating/deactivating  current 
sources.)  Therefore,  transient  thermal  effects  need  to  be 
explored  on  analog  circuit  stand  point,  and  in  order  to  do 
that,  it  is  important  to  have  the  simulator  available. 

Thermal  effects  are  far  more  significant  in  SOI, 
dielectrically  isolated,  and  trench  isolated  transistors. 
Devices  built  with  these  processes  tend  to  be  used  in  many 
high  performance  analog  circuits  because  of  the  advantages 
provided  by  them  (e.g.,  lower  collector- junction  capacitance 
and  lower  isolation  leakage  currents  at  high  temperature) . 
Measurements  indicate  that,  for  dielectrically  isolated 
devices,  thermal  effects  are  about  four  times  more 
significant  but  disappear  at  lower  frequencies  than  same- 
dimension  junction-isolated  devices.  Modeling  thermal 
impedance  for  these  devices  will  be  important  and  very  useful 
for  both  analog  and  digital  circuit  simulations. 
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